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1.  INTRODUCTION 


Background 

Tissue  blood  perfusion  is  a  fundamental  measurement  in 
physiolgy  that  affects  the  entire  spectrum  of  medical  practice 
and  research  [1] .  Many  different  techniques  have  been 
developed  to  measure  tissue  blood  flow.  Two  common  methods 
involve  indicator  dilution  and  radio-labeled  microspheres  [2] . 
These  measurements  are  complicated  and  cannot  be  repeated  at 
frequent  intervals. 

A  different  category  of  blood  perfusion  measurements 
involves  the  use  of  thermal  techniques.  These  methods  have 
the  potential  to  overcome  the  limitations  of  the  indicator 
dilution  and  microsphere  techniques.  However,  roost  thermal 
methods  have  the  drawback  of  comparing  the  heat  dissipation 
against  a  value  of  thermal  conductivity  for  nonperfused 
tissue.  This  requires  either  a  suspension  of  the  blood  flow 
to  the  tissue  or  the  use  of  tabulated  values  for  the  type  of 
tissue  being  e>:amined.  Suspension  of  blood  flow  is  traumatic 
to  tissue  and  might  mean  sacrificing  the  animal,  while  using 
tabulated  values  can  lead  to  inaccuracies  due  to  assumptions 
made  in  creating  the  table.  In  addition,  the  theoretical 
basis  for  some  of  these  methods  is  open  to  question,  since  the 
volume  of  tissue  being  sampled  is  not  much  larger  than  the 
probe  itself,  and  probe  size  and  shape  have  been  shown  to 
alter  results  [ 1] . 


The  thermal  pulse-decay  method  being  developed  at  the 
University  of  Illinois  by  Dr.  Kenneth  R.  Holmes  and  Dr. 
Michael  M.  Chen  is  a  thermal  method  that  overcomes  all  of  the 
problems  described  above.  This  method  consists  of  inserting  a 
small  thermistor  into  the  tissue  of  interest.  A  known 
quantity  of  heat  is  deposited  in  the  tissue  when  current  flows 
through  the  thermistor.  The  thermistor  is  then  used  to 
measure  post-pulse  tissue  cooling,  from  which  local  perfusion 
and  thermal  conductivity  can  be  determined  [1]. 

The  basis  of  these  calculations  arises  from  the  heat 
dissipation  mechanisms  working  in  the  tissue.  With  no  blood 
perfusion,  the  primary  decay  in  temperature  is  due  to  thermal 
diffusion.  This  temperature  decay  can  be  mathematically 
modeled  as  a  decaying  power  series  [1] .  Blood  perfusion  is  a 
thermal  transfer  mechanism  that  can  be  modeled  as  an 
exponential  decay  of  temperature.  ^  vivo  tissue  cooling 
typically  includes  both  of  these  mechanisms  and  can  be  modeled 
as  a  product  of  the  exponential  and  power  series.  From  the 
shape  of  this  curve,  the  thermal  conductivity  and  the  local 
blood  perfusion  can  be  determined  [1]. 

The  advantages  of  this  method  are:  it  provides  an  absolute 
measurement  the  volumetric  perfusion  rate  (ml  blood/ml 

tissue  sec)  without  requiring  calibrations  or  stop-flow 
measurements;  the  sampling  volume  is  considerably  larger  than 
the  volume  of  tissue  traumatized  by  the  microprobe;  the  probe 
shape  and  size  do  not  affect  the  results;  the  electronics  and 
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calculations  are  extremely  simple;  the  increase  in  tissue 
temperature  is  usually  only  0.5  C;  the  small  diameter  of  the 
probe  causes  minimal  trauma  in  the  tissue  under  examination 
[3,4]  . 

Some  drawbacks  of  this  method  are  that  it  does  require 
insertion  of  a  probe,  it  yields  point  information,  and  it  is 
an  indirect  measurement  process. 

Problem  and  Scope 

Small  thermistor  beads  are  fabricated  into  needle-like 
probes  to  aid  in  insertion  into  the  tissue  and  minimization  of 
trauma  in  the  tissue.  The  thermistor  is  incorporated  into  a 
bridge  circuit  that  performs  both  heating  and  measurement 
roles.  The  bridge  output  is  used  to  generate  cooling  curves 
on  a  strip  chart  recorder.  Computer  programs  have  been 
written  to  analyse  these  cooling  curves  for  thermal 
conductivity  and  perfusion  rate. 

The  problem  with  this  procedure  is  that  the  data  must  be 
read  point  by  point  from  the  the  cooling  curves  and  typed  into 
the  computer.  This  process  is  both  time  consuming  and 
inaccurate.  Additionally,  since  the  data  is  analyzed  after  the 
experiment  has  been  completed,  the  operator  has  no  opportunity 
to  modify  the  experimental  parameters. 

The  solution  to  these  problems  is  to  automate  the  control 
of  the  experiment,  the  acquisition  and  storage  of  the  data,  and 
calculations  performed  on  the  data.  This  report  describes  a 
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computer  based  hardware  and  software  system  that  will  not  only 
solve  the  problems  of  entering  cooling  curves  into  the 
computer  and  freeing  the  operator  from  continuous  supervision 
of  the  system,  but  should  also  give  the  operator  the  advantage 
of  having  real-time  blood  flow  data  to  tailor  each  succeeding 
measurement.  Figure  1  is  a  block  diagram  of  the  proposed 
system. 

Presentation 

Chapter  II  describes  the  characteristics  of  the  thermistor 
bridge  circuit,  the  governing  equations,  and  the  desired 
specifications  for  the  proper  operation  of  the  control  system. 

Chapter  III  describes  the  hardware  used  by  the  system 
while  Chapter  IV  documents  the  system  software. 

Chapter  V  presents  the  results  of  experimental  testing  of 
the  entire  system. 

Chapter  VI  contains  the  conclusions  about  the  present 
system  and  recommendations  for  future  systems. 

Software  flowcharts  and  a  complete  program  listing  are 
contained  in  the  appendices. 
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II.  THERMAL  PULSE-DECAY  METHOD 


Descr iption 

The  bridge  circuit  containing  the  thermistor  probe  is  a 
very  simple  balancing  circuit  with  one  modification  that 
allows  current  to  pass  through  the  probe  during  heating.  This 
circuit  is  shown  in  Figure  2.  The  balance  resistor,  is 

used  to  adjust  the  bridge  output.  A  five  Volt  signal 
impressed  at  can  drive  the  transistor  switch  for  the  heater 
relay.  The  duration  of  this  signal  determines  the  length  of 
the  heating  pulse  and  thus  the  energy  delivered  to  the  tissue. 

Simple  circuit  analysis  of  the  bridge  yields  the  following 
relationship  between  the  output  voltage,  E,  and  the  probe 
resistance,  Rp. 

V 

Then,  once  the  resistance  of  the  thermistor  is  known,  the 
temperature,  T,  can  be  found  as: 

T  =  A  -  B  In(Rp)  (2) 

where  A  and  B  are  calibrated  constants  of  the  particular 
thermistor  probe  in  use  [1] . 

The  temperature  rise  due  to  tissue  heating  causes  a 
corresponding  increase  in  the  output  voltage.  Thus,  even 
though  a  nonlinear  relationship  exists  between  voltage  and 
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Figure  2.  Thermistor  Bridge  Circuit 
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temperature,  a  cooling  curve  obtained  by  plotting  the  ouput 
voltage  of  the  bridge  looks  similar  to  a  cooling  curve  that 
plots  actual  temperature. 

For  a  cooling  curve  expressed  in  Volts,  the  typical 
experiment  will  have  about  a  seven  millivolt  range.  In 
addition  to  the  actual  cooling  part  of  the  curve,  important 
information  is  contained  in  the  temperature  values  immediately 
prior  to  heating.  This  pre-sample  period  can  be  used  to 
project  a  baseline  for  use  in  normalizing  the  cooling  part  of 
the  curve  to  temperature  drifts.  A  sample  cooling  curve  is 
shown  in  Figure  3. 

System  Specifications 

The  following  is  a  list  of  desired  specifications  for  the 
automated  control  and  data  acquisition  svstem. 

1.  Deliver  a  square  wave  five  Volt  pulse  to  each  of 
six  bridge  circuits  to  control  the  heating  cycle. 

2.  Vary  the  duration  of  each  pulse  length 
independently  over  a  range  of  zero  to  twenty  seconds 
in  one-tenth  second  intervals. 

3.  The  duration  time  of  each  heat  pulse  must  ba 
stored  for  computation  purposes. 

4.  Monitor  and  record  voltage  output  of  six 
thermistor  bridges  with  a  resolution  of  three 
millivolts  and  a  range  of  minus  five  to  five  volts. 

5.  Allow  for  sampling  rates  up  to  sixty  Hertz  in 
each  of  six  channels. 


6.  Channel  selection  may  be  varied  without  loss  of 
information.  Channels  need  not  be  operated  in 
sequential  order. 

7.  Stored  data  may  be  recalled  for  display  on  a 
video  terminal  or  printed  on  hard  copy. 

8.  System  must  include  a  file  management  system  for 
stored  data  files. 

9.  Calculations  may  be  performed  while  system  is 
not  being  used  to  measure  perfusion. 

10.  Calculate  and  display  steady  state  body 
temperature  as  calculated  from  thermistor  data. 
Calibrate  system  for  each  thermistor. 

11.  If  input  voltage  exceeds  an  operator  selected 
value,  system  will  issue  an  audible  alarm.  Alarm 
may  be  disabled. 

12.  Repeat  measurements  (delivery  of  pulses) 
automatically. 
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III.  HARDWARE 

The  system  hardware  is  composed  of  three  main  units. 
These  are  the  thermistor  bridge,  the  signal  conditioner,  and 
the  computer.  Figure  4  is  a  block  diagram  of  the  entire 
computer  system,  with  the  computer  separated  into  its  major 
components. 

Bridge  Circuit 

The  bridge  circuit  is  described  in  Chapter  II. 

Signal  Conditioning  , 

Preliminary  tests  of  the  bridge  circuit  showed 
considerable  noise  to  be  present,  which  appeared  in  two  basic 
forms.  The  most  prevalent  noise  was  sixty  cycle  power  line 
interference,  which  had  a  magnitude  up  to  four  millivolts, 
almost  as  large  as  the  desired  signal.  The  second  kind  of 
noise  was  a  very  high  frequency  noise  whose  source  could  not 
be  identified. 

In  addition  to  the  noise,  the  signal  strength  was  on  the 
order  of  millivolts,  too  weak  to  be  used  as  an  input  signal  to 
the  analog  to  digital  converter,  which  measures  in  volts.  The 
signal  conditioning  block,  shown  in  Figure  5,  is  used  to 
correct  for  these  two  problems.  First,  the  signal  was 
amplified  by  a  differential  laboratory  amplifier  with  a  gain 
of  1000.  The  input  signal  for  a  typical  experiment  has  a 
range  of  about  seven  millivolts,  with  both  positive  and 


Figure  4.  System  Block  Diagram 
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The  Parallel  Input-Output  unit  is  a  DEC  DRVll.  It 
consists  of  16  separate  input  and  output  lines,  along  with 
appropriate  control  lines.  In  this  application  the  unit  is 
used  to  control  the  heat  pulses  and  light  appropriate 
indicator  lights  on  the  bridge  unit.  In  future  versions,  it 
can  be  used  to  monitor  switch  positions  on  the  bridge  unit. 
The  pin  assignments  for  this  board  are  contained  in  Appendix 
A.  More  information  on  the  operation  of  this  board  is 
contained  in  Reference  [6]. 

The  Analog  to  Digital  converter  is  an  ADAC  1030  which  has 
eight  differential  inputs  each  with  a  range  of  minus  five  to 
plus  five  volts  and  a  programmable  gain  of  1,  2,  5,  or  10.  In 
future  versions  of  the  bridge,  this  gain  can  be  used  as  part 
of  the  1000  gain  of  the  signal  conditioner.  With  the  full 
gain  of  1000,  the  ADC's  12  bit  resolution  is  able  to  look  at 
signals  from  minus  five  millivolts  to  five  millivolts  in  2.5 
microvolt  increments.  This  meets  the  original  specification 
for  resolution.  Timing  for  the  sampling  of  the  ADC  is 
accomplished  with  the  real  time  clock  of  the  LSI  11  computer. 
This  allows  sampling  rates  as  high  as  sixty  hertz  for  each  of 
six  channels.  Additional  inforraat’on  on  the  operation  of  the 
ADC  can  be  found  in  Reference  [7] . 

The  Central  Processing  Unit  is  a  DEC  LSI  11/23  CPU  with 
Memory  Management  Unit  (MMU) ,  model  KDFll-AA,  along  with  the 
floating  point  hardware  option,  model  KEFll-AA.  The  MMU  will 
allow  full  and  efficient  use  of  the  Chrislin  memory  board. 
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The  speed  and  power  of  the  LSI  11/23  CPU  with  floating  point 
should  allow  close  to  real  time  calculaton  of  the  desired 
parameters  after  the  data  has  been  collected. 

The  Chrislin  Cl-1123  memory  provides  a  full  256  kilobytes 
of  random  access  memory.  In  addition  to  containing  the 
operating  program,  it  allows  for  storage  of  quite  large  data 
arrays  during  sampling  that  can  be  transferred  later  to  disk. 

The  Data  Systems  Design  470  disk  drive  provides  two 
megabytes  of  on  line  storage.  Drive  one  will  be  used  as  a 
system  disk  containing  the  source  program  and  monitor,  while 
drive  two  will  be  used  to  interchange  data  disks. 

The  serial  interface  is  a  DEC  DLVll-J,  which  has  four 
independent  communication  channels.  One  of  these  is  used  for 
the  system  terminal.  The  others  could  be  used  for  a  modem,  a 
remote  terminal,  or  a  line  printer. 

The  system  terminal  is  a  Lear  Siglar  ADM5  with  a  512 
retrographics  board.  This  provides  graphics  capability  that 
is  software  compatible  to  the  Techtronics  4010  terminal.  The 
main  use  of  the  graphics  will  be  to  display  various  cooling 
curves  for  visual  inspection.  Additional  information  is 
contained  in  References  [8]  and  [9]. 

The  GPlOO  graphics  plotter  is  connected  directly  to  the 
ADM5  terminal.  It  provides  hardcopy  for  any  of  the  graphics 
sent  to  the  terminal  or  functions  as  a  line  printer  for 
listing  programs  or  data  sets.  It  is  further  described  in 
Reference  [10] . 


I 
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IV.  SOFTWARE 

The  system  software  is  written  in  FORTRAN  and  resides  on 
the  system  disk.  Both  the  user  and  the  operating  program  have 
access  to  an  extensive  library  of  programs  written  for  the  DEC 
LSI  11  computer  system.  These  include  the  RTll  monitor,  a 
text  editor,  a  disk  handler,  linkers,  compilers,  various 

I  input-output  routines,  and  a  system  library  of  FORTRAN 

callable  subroutines.  A  full  description  of  these  programs  is 

f 

•  contained  in  Reference  [11]. 

The  FORTRAN  operating  program  contained  on  the  system  disk 
allows  the  user  to  interact  with  the  system  to  configure  the 
experiment  and  data  collection  in  an  easy  and  flexible  manner. 
In  addition,  it  allows  calculation  routines  to  be  performed  on 
the  data  and  several  modes  of  automatic  operation.  The  user 
also  has  the  ability  to  list  or  graph  any  data  file. 

This  chapter  contains  a  complete  description  of  the 
operating  program.  Figure  6  is  a  simplified  flowchart  of  the 
main  program.  Flowcharts  of  the  subroutines  and  a  complete 
listing  of  the  FORTRAN  code  are  contained  in  Appendix  C. 

Main  Program 

The  program  begins  with  an  initialization  routine  that 
sets  each  variable  to  a  valid  starting  value.  As  the  program 
■  is  used,  the  initial  values  can  be  modified  to  reflect  more 
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The  terminal  is  then  filled  with  the  current  values  of  the 
program  variables  and  a  list  of  the  commands  available  to  the 
user.  An  example  of  this  is  shown  in  Figure  7,  The  commands 
ate  single  letters  that  correspond  to  the  function  or 
parameter  to  its  right.  Each  second,  the  current  time  is 
displayed,  all  six  bridge  outputs  are  sampled  and  displayed  in 
both  millivolts  and  degrees  Celsius,  and  the  time  until  the 
next  automatic  operation  is  updated  and  displayed.  While  in 
this  mode,  the  software  continuously  checks  for  user  input  or 
the  conditions  necessary  to  cause  an  automat..c  operation. 

Each  user  input  causes  the  main  program  to  call  an 
appropriate  subroutine  which  actually  performs  the  command.  A 
summary  of  each  of  the  available  commands  is  shown  in  Table  I. 
The  following  is  a  description  of  each  command  along  with  its 
use  and  restrictions. 

Command  '0' 

This  command  allows  the  user  to  choose  which  channels  will 
be  active  when  either  SAMPLE  or  CALCULATE  is  performed.  Each 
channel  can  be  independently  turned  on  or  off  with  no 
restrictions  as  to  which  channels  must  be  used  for  a  given 
number  of  probes.  The  software  for  command  '0'  is  contained 
in  subroutine  UPDATE. 

Command 

This  command  allows  the  filename  of  any  channel  to  be 
changed.  Each  filename  must  be  constructed  of  three 


20 


TIME  13-  3'  1«  0«Te  S^lE^Sa  DATA  VERSION  ONE 


CHANNEL 

sample  status 

CALCULATE  STATUS 

PTLEHAHE 

PROBE 

CUHRCHT  millivolts 
CURRENT  TEMPERATURE 
LAST  TEMP 

HEAT  <  sec  ■  TICffS  ) 


PRESAMPLE 

SAMPLE 
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ivqlts  b.sm  a.iaa  a.S44  i  1  lit  I  eae 

-PRATURE  34  43Z  35.585  5.888  8  888  «  5  22S 

ERATURE  34.^  8.888  8^888  8^888  8^888 

DURATION  <seC'TIC>  FREQUENCY  CHE)  REAOINCS  PER^O 

8  ««  if®  \ 

15.  a  88  988  1 


L  UPPER  ALARM  LIMIT  f  MV ) 
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e  EXPERIMEHT  INTERVAL  CMIH-SEO  8  • -*5 
EXPERIMENT  HEPETITIOMS  * 

NEXT  EXPERIMENT  < MIN -SEC)  8*  8 

PLEASE  ENTER  YOUR  COMMAND 


AUTOMATIC  RUN 

BREAK  FROM  AUT0 

CALCULATE 

LOOK  AT  DATA  FILE 

SAMPLE 

QUIT 


Figure  7.  Video  Terminal  Presentation  of  Program  Status  and 


Corranands 
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Table  1.  Summary  of  Commands 

Command  Description 

0  Change  Channel  Activity  Status 

F  Change  Data  Filename 

P  Change  Probe,  Bridge  Balance,  or  Descriptive 

Text 

H  Change  Heat  Pulse  Duration 

R  Change  Duration  and  Frequency  of  Sampling 

L  Change  Alarm  Limits 

E  Change  Experiment  Repetition  Parameters 

A  Enter  Automatic  Operation 

B  Break  Prom  Automatic  Operation 

C  Perform  Calculations 

D  Display  Data  Files 

S  Perform  Heat  Pulse  and  Sampling 

Z  Exit  Program 
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alphanumeric  characters  followed  by  three  numeric  characters. 
This  allows  the  user  to  describe  the  experiment  using  three 
letters  and  then  determine  the  experiment  repetition  using  the 
three  numerals  remaining.  This  filename  is  used  by  both 
SAMPLE  and  CALCULATE.  Subroutine  UPDATE  performs  the  command 
'S'. 


Command  'S' 

This  command  allows  the  user  to  change  the  probe  number, 
probe  calibration  data,  bridge  balance  conditions,  and 
descriptive  text  of  any  channel.  When  a  probe  is  changed,  the 
software  reads  in  the  calibration  data  for  the  new  probe  from 
a  disk  file  named  PROBE.DAT.  This  file  is  maintained  by  the 
auxiliary  program  PROBEC.POR  which  is  described  later.  The 
calibration  data  is  part  of  the  data  file  stored  when  a  sample 
is  performed.  It  is  needed  to  calculate  absolute  temperature. 

Changing  the  bridge  balance  conditions  requires  the  user 
to  enter  the  new  resistance  of  the  potentiometer.  This  value 
is  also  needed  for  absolute  temperature  calculation,  and  is 
part  of  the  data  file  stored  by  SAMPLE. 

The  descriptive  text  can  be  used  to  record  probe 
placement,  experiment  objective  or  other  pertinent 
information.  It  is  also  part  of  the  data  file  stored  by 
SAMPLE . 

Subroutine  UPDATE  performs  the  command  'S'. 


23 


Command  'H' 

This  command  allows  the  uset  to  modify  the  duration  of  the 
heat  pulse  that  is  applied  to  the  probe.  Each  channel  can  be 
varied  independently  in  increments  of  one-sixtieth  of  a 
second.  A  duration  of  zero  is  allowed  and  simply  means  that 
that  channel  will  not  be  pulsed  during  sampling.  This  allows 
probes  to  be  used  to  determine  heat  patterns  delivered  by 
other  probes.  Subroutine  UPDATE  performs  the  command  'H'. 
The  heat  pulse  duration  is  part  of  the  data  file  stored  by 
SAMPLE. 

Command 

This  command  allows  the  user  to  define  the  sampling  rate 
and  duration  for  an  experiment.  The  samples  taken  before  the 
heat  pulse,  the  presample  period,  can  be  configured 
independently  of  the  samples  taken  after  the  heat  pulse.  The 
user  inputs  the  desired  time  duration  and  frequency  of 
sampling.  The  program  calculates  and  displays  the  actual 
number  of  readings  to  be  taken  and  the  time  period  (in 
one-sixtieth  of  seconds)  between  samples.  These  values  are 
part  of  the  data  file  stored  by  SAMPLE.  Subroutine  UPDATE 
performs  the  command  '“R"’. 

Command  'L' 

Thisoommand  lets  the  user  set  limits  for  an  audible  alarm 
that  monitors  the  input  voltages  of  the  channels  that  are 
active  for  sampling.  If  the  input  voltage  is  not  within  the 
specified  range,  the  program  sounds  the  bell  on  the  terminal 
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each  second.  With  the  correct  limits,  the  alarm  can  be  used 
to  notify  the  user  when  a  bridge  needs  rebalancing  before 
starting  an  experiment.  Subroutine  UPDATE  performs  the 
command  'L". 

Command  'E' 

This  command  allows  the  user  to  define  the  parameters  for 
operating  under  automatic  control.  The  time  interval  between 
experiments,  the  number  of  experiments  to  be  performed  and  the 
time  until  the  first  experiment  can  all  be  set  to  any  value. 
Subroutine  UPDATE  performs  the  command  'E' . 

Command  'k' 

This  command  enables  the  automatic  operation  of  the 
program,  as  defined  by  the  automatic  control  parameters.  It 
is  performed  by  the  main  program  and  subroutine  AUTO.  When 
called,  subroutine  AUTO  will  call  SAMPLE  and  CALCULATE  as 
needed,  increment  the  filenames  of  the  active  channels  and 
update  the  automatic  operation  parameters.  Using  the  'A' 
command,  the  operator  can  free  himself  from  having  to 
continuously  monitor  and  initiate  experiments  or 
calculations . 

Command  'B' 

This  command  disables  automatic  operation  of  the  program. 
The  time  until  the  next  experiment  will  continue  to  count 
down,  but  AUTO  will  not  be  called  if  the  count  reaches  zero. 
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Using  the  'B'  command,  the  user  can  exit  automatic  control, 
change  a  system  parameter,  reenter  automatic  control  and  not 
lose  the  correct  spacing  between  experiments.  The  main 
program  performs  the  'B'  command. 

Command  'C' 

This  command  performs  the  desired  calculation  on  each 
channel  that  is  active  for  calculate.  Subroutine  CALC 
performs  this  command  but  does  not  perform  the  calculations. 
CALC  reads  and  prepares  the  data  file  for  a  subroutine  called 
CRUNCH,  which  is  to  perform  the  actual  calculations.  CRUNCH 
can  either  be  a  thermal  conductivity  or  blood  perfusion 
routine  that  can  be  linked  to  CALC  for  use. 

Command  'D' 

This  command  allows  the  user  to  examine  a  data  file  by 
listing  or  graphing.  Listing  a  data  file  consists  of 
labeling,  formatting  and  printing  the  system  configuration  at 
the  time  of  sampling  and  the  entire  presample  and  sample  data. 

The  graphics  package  plots  the  cooling  curve  and  labels  it 
with  pertinent  information  from  the  data  file.  The  user  has 
the  options  of  changing  the  size  and  position  of  the  graph  as 
well  as  expanding  the  time  scale  to  display  only  a  portion  of 
the  cooling  curve. 

Subroutine  DATA  performs  the  command  'D'. 


Conunand  'S 


This  conunand  causes  the  computer  to  perform  the  sampling 
as  defined  by  the  system  parameters.  For  each  active  channel 
the  appropriate  presample  is  collected,  the  probe  is  heated, 
the  desired  sample  data  is  collected,  and  a  data  file  is 
created  on  disk.  The  format  of  the  data  file  is  in  appendix 
B.  It  consists  of  all  possible  information  that  might  be 
desired  at  a  later  time  to  analyse  the  results  of  the 
experiment.  Subroutine  SAMP  performs  the  command  'S'. 

Command  '2' 

This  command  terminates  the  operation  of  the  program  and 
returns  control  to  the  RTll  monitor. 

Auxiliary  Program 

The  auxiliary  program  PROBEC.FOR  is  used  to  maintain  the 
data  file  that  contains  the  current  probe  calibration 
constants.  The  user  has  the  option  of  entering  calibration 
data  for  a  new  probe  or  changing  calibration  data  for  an 
existing  probe  in  the  file.  Each  time  new  data  is  entered 
into  this  file,  the  user  can  also  enter  the  date  of 
calibration.  The  probe  calibration  and  the  date  are  stored  in 
this  for  use  by  the  system  software.  Specifications  for  the 
data  file  PROBE.DAT  are  contained  in  appendix  B. 
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V.  RESULTS 


The  data  aquisition  system  has  been  successfully  used  in  a 
variety  of  configurations. 

Figures  8  through  14  are  actual  output  from  the  system. 
Figure  8  shows  the  cooling  curve  of  a  probe  that  was  placed  in 
a  100  ml  solution  of  glycerin  at  38  degrees  Celsius.  The 
presample  data  is  shown  from  minus  five  to  minus  two  seconds. 
The  plot  disappears  from  minus  two  to  zero  seconds,  since  that 
is  during  the  heat  pulse  and  the  bridge  output  is  not  valid. 
The  sample  data  is  shown  from  zero  to  fifteen  seconds. 

This  curve  shows  that  the  probe,  the  bridge,  the  signal 
conditioning,  and  the  analog  to  digital  converter  all  worked 
properly  to  deliver  the  appropriate  data  to  the  computer.  It 
also  shows  that  the  software  is  capable  of  procuring  the  data, 
storing  it,  and  retrieving  it  for  display. 

Figure  9  is  a  cooling  curve  from  a  similar  experiment  but 
the  heating  pulse  has  been  reduced  to  only  one-sixth  of  a 
second.  The  curve  decays  much  more  rapidly  since  less  heat 
was  deposited  by  the  probe.  In  order  to  display  this  data  in 
a  more  visible  manner,  the  software  allows  for  expansion  of 
the  time  scale  so  that  any  portion  of  the  curve  can  be  shown. 
An  example  of  this  capability  is  shown  in  Figure  10  using  the 
same  data  file  as  Figure  9. 

Figure  11  is  from  a  data  file  that  was  created  from  an 
experiment  using  a  live  rabbit.  The  blood  perfusion  in  the 
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kidney  could  be  determined  from  this  data. 

Figures  12  and  13  are  from  data  files  created  from 
experiments  conducted  on  the  liver  of  a  live  rabbit.  The 
probe  for  Figure  12  was  placed  very  close  to  a  large  blood 
vessel  in  the  liver.  This  results  in  a  cooling  curve  that 
decays  quite  rapidly.  Also  present  on  this  curve  is  some 
respiratory  artifact  that  causes  the  curve  to  rise  and  fall 
periodically.  Figure  13  is  from  a  test  on  the  same  liver  but 
the  probe  has  been  moved  away  from  any  blood  vessels. 
Perfusion  is  still  present,  but  at  a  much  lower  rate  and  the 
respiratory  artifact  is  no  longer  visible.  In  order  to 
compare  these  two  curves.  Figure  14  was  generated  by  expanding 
and  overlaying  the  data  from  the  two  files  on  the  same  plot. 
The  difference  between  the  high  perfusion  curve  and  the  low 
perfusion  curve  is  now  seen  as  different  rates  of  decay  in  the 
cooling  curves.  The  higher  perfusion  near  the  large  blood 
vessel  rapidly  removes  heat  from  the  vicinity  of  the  probe, 
resulting  in  a  faster  decay  than  is  seen  in  the  low  perfusion 
case. 
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Figure  13.  Rabbit  Liver,  Low  Blood  Perfusion 
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RABBIT  LIUER  PERFUSION  TEST 
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14.  Rabbit  Liver,  Comparison  of  High  and  Low 
Perfusion 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

An  automated  control  and  data  aquisition  system  for  the 
thermal  pulse-decay  method  of  blood  per fusion  measurement  has 
been  designed  and  tested.  The  system  meets  or  exceeds  the 
desired  specifications,  providing  a  useable  tool  for  the 
continuation  of  research  in  local  blood  perfusion  measurement. 

Recommendations 

Based  on  system  performance  until  this  point,  the 
following  recommendations  are  proposed  as  future  modifications 
to  the  system.  Additional  desired  modifications  will  become 
apparent  as  the  system  is  utilized  more  frequently. 

1)  Incorporate  into  the  auxiliary  program,  PROBEC.FOR, 
the  capability  to  automatically  calibrate  the  thermistor 
probes. 

2)  Enhance  the  data  presentation,  particularly  the 

ability  to  show  more  than  one  cooling  curve  on  the  same  plot. 

3)  Provide  more  disk  file  security  to  prevent  accidental 
overwriting  of  data. 

4)  Provide  for  easy  transfer  of  data  files  to  or  from 
this  computer  and  another. 

5)  Incorporate  computer  routines  to  analyse  the  cooling 
curves  for  be  :h  thermal  conductivity  and  blood  perfusion. 
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APPENDIX  A 

SYSTEM  HARDWARE  CONNECTIONS 


This  appendix  contains  the  wiring  list  for  connecting 
the  DEC  DRVll  parallel  input-output  board  to  the  bridge 
circuit.  Two  twenty- five  wire  ribbon  cables  carry  the 
signals  listed  below  in  Table  2. 


Table  2.  Pin  Assignments  for  Bridge  Control 


Signal 

Channel 

Connector , 

Pin 

Presample  LED 

1 

Jl, 

11 

Presample  LED 

2 

Jl, 

12 

Presample  LED 

3 

Jl, 

13 

Presample  LED 

4 

Jl, 

14 

Prescunple  LED 

5 

Jl, 

15 

Pres^unple  LED 

6 

Jl, 

16 

Heat  Pulse  and 

LED 

1 

J2, 

11 

Heat  Pulse  and 

LED 

2 

J2, 

12 

Heat  Pulse  and 

LED 

3 

J2, 

13 

Heat  Pulse  and 

LED 

4 

J2, 

14 

Heat  Pulse  and 

LED 

5 

J2, 

15 

Heat  Pulse  and 

LED 

6 

J2, 

16 

Saunple  LED 

ALL 

Jl, 

17 
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APPENDIX  B 
DATA  FILE  STRUCTURES 

This  appendix  contains  the  structure  of  the  two  types 
of  data  files  used  by  the  operating  software.  PROBE.DAT 
contains  the  calibration  data  for  the  thermistor  probes  and 
is  outlined  in  Table  3.  The  structure  of  the  data  files 
created  by  SAMPLE  is  outlined  in  Table  4 . 


Table  3.  PROBE.DAT  Structure 


Record 

Variable  Description 

Disk 

Format 

First 

Number  of  Probes 

1*2 

Second 

Probe  Number 

1*2 

Month 

1*2 

Day 

1*2 

Year 

1*2 

A 

R*4 

B 

R*4 

RBO 

R*4 

BETA 

R*4 

Repeating  for  each  probe  in 

the  data 

file . 

Table  4. 

Data  File  Structure 

Variable  Description 

Variable  Name 

Disk  Format 

Data  Version 

1 

1*2 

Time 

ISTIM 

1*4 

Date 

KDATE 

3  1*2 

Channel 

I 

1*2 

Text 

ITEXT 

20  1*2 

Active  Files 

IFILE(l) 

3  1*2 

Active  Files 

IPILE(2) 

3  1*2 

Active  Files 

IFILEO) 

3  1*2 

Active  Files 

IFILE(4) 

3  1*2 

Active  Files 

IFILE(5) 

3  1*2 

Active  Files 

IFILE(6) 

3  1*2 

Probe  Number 

IPROBE 

1*2 

RBO 

PPRBO 

R*4 

BETA 

PPBETA 

R*4 

A 

PPA 

R*4 

B 

PPB 

R*4 

Probe  Calibration  Date 

ICDATE 

3  1*2 

Bridge  Balance  Resistance 

IPPRBB 

1*2 

Bridge  Balance  Time 

IPPTBB 

1*4 

Bridge  Balance  Voltage 

PPVBB 

R*4 

Heat  Pulse  Duration  (Ticks) 

IHEAT 

1*2 

presample  Readings 

IPSR 

1*2 

Presample  Period  (Ticks) 

IPS? 

1*2 

Sample  Readings 

ISR 

1*2 

Sample  Period  (Ticks) 

ISP 

1*2 

Presample  Data 

IDATA 

(IPSR)  1*2 

Sample  Data 

IDATA 

(ISR)  1*2 
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APPENDIX  C 
OPERATING  PROGRAM 

Figure  6  in  Chapter  IV  is  a  flowchart  of  the  main 
program.  This  appendix  contains  a  listing  of  the  FORTRAN 
code  of  the  main  program  followed  by  flowcharts  and  code 
listings  of  the  subroutines - 

Figure  15  shows  the  flowchart  of  subroutine  INIT  (p.  46) , 
Figure  16  shows  the  flowchart  of  subroutine  UPDATE  (p. 49-50), 
Figure  17  shows  the  flowchart  of  subroutine  SAMP  (p.  57) , 

Figure  18  shows  the  flowchart  of  subroutine  CAI.C  (p.  63)  , 

Figure  19  shows  the  flowchart  of  subroutine  AUTO  (p.  65), 

Figure  20  shows  the  flowchart  of  subroutine  DATA  (p.  63) , 

Figure  21  shows  the  flowchart  of  PROBEC.FOR  (p.  74) ,  and 

Figure  22  shows  the  flowchart  of  subroutine  GET  (77) . 
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THIS  IS  THE  MOIN  PROGRAM  FOR  THERMAL  PULSE 
C  STUDIES 

-  IT  UAS  WRITTEN  BY  KURT  L  BAUM 

C  tt  4t  *  COMMON  variables  *  *  :« 

LOCtCALtt  HAZ< 3 HAZA< S > 

INTEGER  IS<23£> 

INTEGER*4  I4<12> 

REAL*4  R4<  fi4  > 

integer  tOATA<iaaQ> 

COMMON  yHAZEL-'HAZ  .  HA2A 
COMMON  ✓VAHBLS^ia.  r4,R-t 
COMMON  -'TEMPS-' I  DATA 
C  t  C  «  ENO  COMMON  BLOCK  XXX 

G  X  X  «  BEGIN  COMMON  OESIGNATICNS  XXX 

INTEGER  IACTVS<  Si.  lACTVCC  S)-  IPROBEC  S  >-  IHEAT' 6  >-  IPPR29C  S > 

INTEGER  IFILE<  S.  3  )  -  ITE.XTf  S-  20  i.  ICDATEC  S  -  3 
INTECERX-4  tPPTBBCS) 

REALX4  PPRaa<  S ) >  PPBETAC  S i - PPA<  S  > - PPSC  6 i - PPVeSC  S i 
REALX4  CTEMP<«J,RLTe«P<E> 

equivalence  < icon, I2<  I  >  ). 

EQUIVALEHCe  C IPROBCC 1 > - r2< S > > - 

EQUIVALEHCe  <  IPPRSBCl  > .  ISC  29 > J . 

EQUrVALEHCe  C IPSP- I2C27> >. 

cauivALEHce  <■  ISP.  rac29  >  >- 

equivalence  <  rCOATEC I  -  I  >, I2<  1 50  >  > - 

EQUIVALENCE  < IFILEC I  -  1  V.  t2< 174  >  >- 

EQUIVALENCE  < lEXREP . I 2< 193 > > 

EQUIVALENCE  <  IPPT8e<  1  > -  1 4<  1  > > 

equivalence  < PPRBa< I >, R4< 1 > >. 

EQUIVALENCE  < PPA< I > . R4< 13 > > - 

EQUIVALENCE  < PPVB8<  1  > - R4< 25  )  > - 

EQUIVALENCE  <  RLTEMPC  I  > .  R4<  37  . 

equivalence  < RLALIM- R4< 44 > >, 

EQUIVALENCE  < AUTOTU > R4< 40 > > . 

C  XXX  ENO  COMON  OESICNATIONS  XXX 

C  XXX  e~GIN  LOCAL  VARIABLES  XXX 
LOGi  ;ALX1  PRNTIM<91 
INTEGER  JTEMP<S) 
integer  XTYPEC30) 

INTEGERX4  KTIME 
INTEGER  KTIM<2) 

REALX4  CV0LT< S  ), TEMP 

EQUIVALENCE  < KT I  ME- KT I M< I  > > 

C  XXX  ENO  LOCAL  VARIABLES  X  X  X 

C  XXX  BEGIN  CODE  SEGMENT  XXX 
HAZ<  2  :>'>''20e 
HA2C  3  >-"200 
HA2A< 1 >-27 
HAZAC  2  >-01 
HAZA<  5  >-"200 
CALL  IN  IT 

PRNTIM< 9 >-*200  i 

C  X  X  X  PRINT  PROGRAM  IDENTIFICATION  X  «  X 

type  30 

HAZC I >-2S  , 

CALL  PRINTCHAZ> 

HAZAC 3 >-49  \ 


<  IACTVSC 1  > .  I2C  2  >  > 

•:  iHEATc  1  :■ .  lac  1 4 
c  iPSR.  i2<  2s  :■  :■ 
c  isR.  I2C  aa  >  > 

C  ITEXTC  1  -  I  >,  I2C33  >  > 
< lACTvcc 1 >, I2C lea  >  > 
C  lAUTOC.  I2C  192  >> 


C  PPBETAC  I  1 , R4C  7  >  > 
C  PPBC  1  >- R4<  19  >  > 

C  CTEMPC  1  >- R4C  31  >  > 
C  RUALIM- R4C  43  >  > 

<  AUTOTN , R4C  45  >  > 

C  EXINT- R4<  47  >  > 


iM  O 


42 


HAZA< 

CALL  PRINT? HAZA> 

TYPE  1 

1  FORMAT?  • THERMAL  PULSE  DECAY  CONTROLLER') 
TYPE  30 

TYPE  2 

2  FORMAT? T1 3.  ' WRITTEN  8Y  KURT  L.  9AUM '  ) 
ACCEPT  3 

3  FORMAT?  A ) 


« 


100 


103 


107 
1  10 


HE 
113 
120 
133 
13S 
137 
133 
1  40 
1-42 


143 


*  K  PRINT  SMORGASBQARO  «  X  « 

TYPE  30 

HAZ? 1 )-2S 

CALL  PRINT? HAZ) 

CALL  lOATE? KMON. KOAY- KYEAR > 

ENCOOE?  3.  100. PRHTIM  >  KMOH. KOAY . KYEAR 
format?  I2,'/’‘,I2,‘^'.I2> 

CALL  PRINT? 'TIME' ) 

HAZA?  3 )>32 
NAZa?  4 )»3S 
CALL  PR  I  NT? HAZA) 

CALL  PRINT? 'DATE' ) 

HAZA?  4  )-37 
CALL  PRINT? HAZA) 

CALL  PRINT?PRNTIM> 

HAZA?  4 )-77 
CALL  PRINT? HAZA) 

CALL  PRINT? 'OATA  VERSION  ONE') 

TYPE  103. 1.2. 3. 4. 3. e 

FORMAT?'  CHANNEL  -.S??)?,!!)) 

□O  107  I-l.E 
JTEMP? I )«IACTVS? I ) 

IF? lACTVS? I ) . NE  a )  GOTO  107 
JTEMP?!)-' 

CONTINUE 
TYPE  110. JTEMP 

FORMAT? 'O  SAMPLE  STATUS  •,S?SX.A2)) 

00  112  1-1.6 
JTEMP?  I )-IACTVC?  I  ) 

IF?  lACTVC? I  )  . NE . 0  )  GOTO  ll2 
JTEMP? I)-' 

CONTINUE 
type  1 13. JTEMP 

FORMAT?'  CALCULATE  STATUS  '  . 6? 6X . A2  )  ) 

TYPE  120.  ?  ?  I  FILE?  I.K>.t?-1.3).I-1.6> 

FORMAT?  'F  filename  .eC2'X.3A2)) 

TYPE  133.IPR0eE 

FORMAT?  'P  PROBE  '.610) 

TYPE  136 

FORMAT?'  CURRENT  MILLIVOLTS') 

TYPE  137 

FORMAT?  ■  CURRENT  TEMPERATURE'  ) 

TYPE  138.RLTEMP 

FORMAT?'  LAST  TEMP  ' . SFB  3  ) 

TYPE  140.?  IHEAT?  I  ).-S0.  IHEAT?  I  )-IHEATC  I  >.'60*60.  I«1 . 6  > 
FORMAT?  'H  HEAT  ? SEC 'TICKS)  '  . 6?  1 3 .  '  '  '  .  1 2  )  ) 

TYPE  142 

FORMAT?'  DURATION  ? SEC  TIC)  FREOUENCY 

1  READINGS  PERIOD') 

KTIM?  1  )-0 

KTIM?  2 )-IPSR*  IPSP 

call  CVTTIM?  KTIME. KHRS. KMIN. Ksec. KTICXS  > 

TYPE  143.  KSEC.  KTICKS  .  60.'  I  PSP  .  IPSR-  I  PSP 
FORMAT? 'R  PRESAMPLE  ' . 12. ' ' ' . 12- ' 

1  I Z  -  '  '.14.'  ',12) 


?  HZ  /  •  , 


.<Tri>i<2  >-:sRt!r£? 

CftLL  CVTT  rM<  itTIME.  KHRS  .  k:mIH,!<SEC.KTICK3> 

TYPE  ISa,  KSEC  .  KTICI'S  -  ISP.  I  SR.  ISP 
1=0  rORMPT<  '  SPMP'-E  •  .  13.  '  ‘  .  13.  ' 

1  13-  •  •  .  13  > 

TYPE  30 

TYPE  ISl.RURLIM 

151  FSRMAT<‘t_  UPPER  PUARM  LIMIT  <  MV  >  '  .  FS  . 

1  T>S. 'R  AUTOMATIC  RUN" > 

TYPE  ISa.RLALIM 

153  FORMAT<  '  LOWER  ALARM  LIMIT  < MV >  • . F5 

1  TPS.'B  BREAK  FROM  AUTO •  > 

KTIM<  1  )-:HT<  EX!NT..S3S3fi  B  > 

•<TIM<  3  >-INT<  SXINT-SSS3fi  .  a;»KTIM<  1  >  > 

CALL  CVTTIM<  KT I  ME . KHRS . KM IN . KSEC . KTICKS  > 

TYPE  133. KMIN. KSEC 

133  FORMAT<  'E  EXPERIMENT  INTERVAL  < MIN  SEC >  '.13. 

1  T43.  'C  CALCULATE*  > 

TYPE  ISa.IEXREP 

1S0  FQRMAT<  •  EXPE.'^IMEHT  REPETITIONS  '  .  I 

1  T4S.  *0  LOOK  AT  DATA  FILE*  > 

TYPE  ISS 

IS3  FORMAT< *  NEXT  EXPERIMENT  < MIN  SEC/  0'  0 

1  T-iS.  *  S  SAMPLE  *  ) 

TYPE  1S« 

ISfi  FORMAT<  T.t3.  *  2  QUIT*) 

HA2A<  3  >-S3 
HAZA<  A  »33 
CALL  PRINTCHAZA) 

TYPE  190 

190  FORMAT<  *  PLEASE  ENTER  YOUR  COMMAND *.» > 

30  FORMATC  > 

tF<  ICOM.NE.es  )  GOTO  30 
TYPE  30 

CALL  PRINTCHAZA) 

TYPE  MO 

M0  FORMAT<  *  UNDER  AUTOMATIC  OPERATION*  > 

C  *  *  *  CHECK  USER  INPUT  *  *  * 

30  CALL  IPOKE<  *’MM.  **010100  OR  .  IPEEK<  '•MM  >  > 

ITEMP-ITTINR<  } 

CALL  IPOKEC  **MM.  **  I .  AND  IPEEKC  "MM  >  > 

IF  <ITEWP.GE.O)  GOTO  300 

C  %  ■*  *  CHECK  FOR  JUMP  TO  AUTO  X  *  * 

I F< r COM . EG. 6S. AHO. AUTOTN. EG. 0 . 0>  GOTO  307 

C  *  *  *  CHECK  FOR  TIME  AND  TEMP  UPDATE  t  *  * 

T7  CALL  CTIMKKTIME) 

CALL  C*/TTIM<  KTIME.  KHRS.  KMIN.  KSEC.  KTICKS  ) 

IF  C KSEC . EG . ISEC  1  GOTO  SO 
ISEC'KSEC 

C  *  *  *  REPRINT  TIMES  AND  TEMPS  »  *  * 

ENCODE*:  a.  80.  PRNTIM  )  KHRS.  KMIN.  KSEC 
30  FORMATC 13. *■*, 13.  *.*. 13 > 

C  *  ;«  4t  GET  CURRENT  TEMPS  *  *  « 

K«1 

DO  90  I >1.6 

CALL  IPOKEC * 17S770. K > 

J-IPEBKC  " 17S773  ) 

IF< J . CT . "3777 )  J-J-"1900O 
CVOLTC  I  J-Jz-ZBa  903 
C  *  «  *  309  90S-"3777zS  .  33  JI  *  « 


u  »•!  n  Cl 
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iF< lACTvsc I >  EQ . a >  GOTa  aa 

IFC  CVOLTC  I  ;■  GE  RLALIM  AMO  C'/Ot_T<  I  >  LE  .  RU  AL  1 11  >  CO’'a  88 

HA2C 1 ^-7 

CAUL  PRINT? HAZ> 

38  RB8>FL0AT<  iPPRse?  I  )  >^a  a-Meaa  a 

VRS-l 197S« . as 
RS-ase7S . a 

TEHP»VRS^<  VRS^C  RBB^RS  >*CVOLT<  I  >^1908  >-9S 
C»  *  *  UHERE  1  197SS  .  aa-V8*a2a73-  VB-5  .-faS  VOLTS  t  Jt  * 

CTEMP<  I  >-PPA<  I  )-PPB<  I  )*ALOC<  RB  > 

K»K*"-4aa 
90  CONTINUE 

C  *  *  *  PRINT  TIME  *  *  * 

TYPE  30 
HAZA<  3  >>32 
HAZA<  ■*  :>3a 
CALL  PRINT? HAZA> 

CALL  PRINT? PRNTIM ) 

C  «  *  t  PRINT  CURRENT  MILLIVOLTS  *  «  * 

TYPE  30 
HAZA<  3  >>38 
HAZA?  T ixSS 
CALL  PRINT? HAZA) 

TYPE  SS.CVOLT 
9S  FORMAT? SFB . 3 ) 

C  *  t  *  PRINT  CURRENT  TEMPS  *  *  » 

TYPE  30 

HAZA? 3 >-39 

CALL  PRINT? HAZA) 

TYPE  97.CTEMP 

97  FORMAT? 6F8 . 3  ) 

C  »  »  *  UPOATE  ANO  PRINT  NEXT  AUTO  TIME  *  :*  * 

IF?  AUTOTN  £0.0.0)  GOTO  99 
IF?  lAUTOC . EQ . 0  )  GOTO  9S 
TIME-«3S3S  .  a*KTIM?  1  >-*-KTIM?a) 

IF? AUTOTU . CT . TIME  )  AUTOTU-AUTOTU-S38T000  0 
AUTOTN- AUTOTN-?  T I ME-AUTOTU ) 

IF? AUTOTN . LT . 0 . a )  AUTOTM-0 . 0 
AUTOTU-TIME 
9€  TYPE  30 

HAZA?  3 )-30 
HAZA?  9 )-60 
CALL  PRINT? HAZA) 

KTIM? 1 )-IHT? A0T0THXSS53S. 0 ) 

KTIM?  a  )-rNT?  AUTOTN-S5530  0»(?TIMC  1  >> 

CALL  CVTTIM?  KTIME.  tCHRS.  KMIN,  KSEC .  KTI Cf?S  > 

TYPE  98.  tCMIM,  KSeC 

98  FORMAT? 12.  la  ) 

C  ■*  X  :*  RESTORE  CURSOR  *  *  X 

99  TYPE  30 
HAZA?  3  )-S3 
HAZA? 9 >-39 

CALL  PRINT?HAZA> 

GOTO  30 

*  *  I  PERFORM  RCCIEVeO  COMMAND  *  *  » 

aa  ICOM-ITSMP 

03  TYPE  30 

07  IF  CtCOM.NE.fiS)  GOTO  aiS 

IF  ? lAUTOC . NE . a )  GOTO  219 


tftUTOC-l 

CALL  CTtM<KTIMe> 

AUTOTU-6S33S  0;«KT  I  M<  1  >+KT  t  M<  a  > 

I F<  AUTOTN  £0.00)  CALL  AUTO 
IF  ■:  ICOM  .  EQ  .  67  >  CALL  CALC 
IF  < icoM  ea  sa>  call  data 
IF  f ICOM . EQ . 90 >  GOTO  330 
IF  (icon  £0.93 >  CALL  SAMP 

IF  <  ICOM , GT . 68 . ANO  ICOM . LT . 33  )  CALL  UPDATE 
GOTO  3 

STOP •  END  OF  program 

END 


Set  All 
Common  . 
Vairiables 
to  Zero 


Figure  15.  Flowchart  of  Subroutine  INIT 
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iUBROOTTNE  IMI” 

TKrs  sceRouTiNE  iMiTiHriiLicss  the  svsreri  poramaters 

«  K  *  common  VARIAELES  t  S  « 

LOCrCftLtl  HR2C 3 >, HftZA< 5 > 

TMTECEH  t2<2SE) 

rnTecc»«-t  r-KiEj 

REai.*-t  R-^Cfi-O 

tNTECSR  rORTRC 1800) 

COMMON  v'HAZEL/'HAZ,  HftZA 
COMMON  ^V#«Rei.E^  12  .  t-».  9-t 
COMMON  ^'TEMPS^  rORTfl 
(  (  «  EHO  COMMON  8(.QCK  Ic  X  « 


X  8ECIN  COMMON  QES  I  ONfiT  I ONS  XXX 

INTEGER  rRCTVSi:  6  ).  I  ACTVCC  S  V  .  IPROBEC  6  '>  .  tHEAT<  S 
INTEGER  tFiL£<8. 3  i.  ITEXTCS, 20 ).  ICOATE: S. 3  > 


IpoftB9<  S  '/ 


eauiVAcEHCE 

eaUIVALEMCE 

EaiJIVACEHCe 

eaurvALEHCE 

equivalence 

EQUIVALENCE 

EQUIVALENCE 

equivalence 

EQUIVALENCE 

equivalence 
equivalence 
equivalence 
X  X  «  end  comoh 


C  ICQM, I2<  t  >  > . 

<  IPROBEC 1 >,  I2<  8  >  >- 

<  IPPRBBC 1  ).  I2<  20  > 

< IPSP, I2<  27 ) >- 

<  ISP,  I2<  23  >  >, 

<  ICQATEC  l,l>.la<l50>>. 
< IFILEC 1, I >, I2C 174 >>- 

<  I EXREP,  I2< 193 >  > 

<  PPVBBC 1  > , R4<  25  >  > , 

<  RLTEMPC  1  ) ,  R4<  37  i , 
<RLALIM,R4<44>>. 

< AUTOTU, R4<4«>  >, 
DESIGNATIONS  XXX 


c  iactvsc  1  .  I2C  2  V  > 

C  I  NEAT-:  1  >,  :2C  14  >  i 

<  ipsr. :2<ze > > 

(  ISR  ,  i2<  23  >  ) 

C  ITEXTC  I  -  1  >  ,  I2C30  i 

<  lACTVCC  1  12'  ’ 

<  I autqc , lac  1 9 

C  CTEMPC  1  > . R4C  31  > 1 
C  R'JAL  IM,  RT<  43  >  > 

C  AUTOTM- R4<  45  >  > 

<  EXINT, R4<  47>  > 


BEGIN  CODE  SEGMENT  X  X  x 
OO  100  I>1,256 
I2< I )-0 
CONTINUE 

OO  110  Iol,64 
R4< I >-0 . 0 
CONTINUE 

□0  120  I«1 , 12 
14< I )-0 

continue 


□0  5  1-1 , S 
If-  ILE<  I  ,  1  ‘  AA  ' 

, PILEC  I , 3  3- • 00 • 
continue 
IFILEC  1 , 2  >- ■ A0 ' 
IFILEC  2,  2  >-  '  90  • 
IFILEC  3, 2  ' CO ■ 

IFILSC  4, 23-'O0’ 
IFILEC  3- 2  3- ‘ E0 ' 
IFILEC  S, 2  3- • F0* 


(*)  o 
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3ijbrcijt:he  upoi^te 

THIS  SUBROUTINE  ftULOUS  THE  USER  TO  CHWNCE 
PARAMETERS  TO  TAIUCR  SYSTEM  COUP  I  OUR AT  I OM 

C  S  S  t  COMMON  VARIABLES  t  *  « 

LOS  I  CAL'*  1  HAZ<  3  ),  HAZAt  5  > 

IHTEOER  ia<2SB> 

INTEGER  It<12> 

REAL*-*  R-*<  64  ) 

INTEGER  inATA<lS00> 

COMMON  -HAZEL^HAZ- HAZA 
COMMON  ^VAR8LEXia> I4>R4 
COMMON  -'TEMPS^  lOATA 
C  t  *  *  END  COMMON  BLOCK  *  »  » 

C  X  -I  t  BEGIN  COMMON  OES I  GNAT  IONS  «  «  « 

INTEGER  lACTVSC  6  >.  I  ACTVCC  6  >.  IPROBEC  S>.  IHEATCS;. ,  IPPRSSC-5> 
INTEGER  IFILE<S-3>.  I TEXT<  6, 20  >-  ICDATEC  6.3) 

INTECER*4  IPPT8B<6) 

ReAL*4  PPRS0<  S ). PPBETA<  6  >. PPAC  6  >. PPB<  6  >- PPVBBC  6  ) 

REAL«4  CTEMP< 6 ) . RLTEMPC 6 > 

EQUIVALENCE  < ICOM. :2< 1 ) ) . 

EQUIVALENCE  <  I PROBE!  I  ) -  I 2< Q >  )  . 

EQUIVALENCE  < I PPRBB! 1 > . I 2< 20 > > . 

equivalence  <  IPSP,  I2<  27  >  •> - 

EQUIVALENCE  < I SP. I2< 29 > > > 

equivalence  <  ICOATE!  1 .  t  >.  I2C  150  >  >. 

equivalence  < I  FILE! 1.1).I2<174>>, 

EQUIVALENCE  <  I EXREP  -  I 2<  1 93  )  > 

EQUIVALENCE  <  I PPTB8<  1  > .  I 4<  1  >  ) 

EQUIVALENCE  < PPRBBC  1  R4<  1  >  ) . 

EQUIVALSNCE  < PPA<  I  ), R4<  13 > >. 

EQUIVALENCE  < PPVBBC 1 > . R4< ZS  )  > . 

equivalence  < RLTEMP! 1 ). R4< 37 > >- 

EQUIVALENCE  C RLAU I M . R4< 44 >  )  . 

EQUIVALENCE  < AUTOTU . R4< 46 > > , 

C  X  *  »  END  COMON  OESIGNATIONS  S  X  « 

C  XXX  BEGIN  LOCAL  VARIABLES  XXX 
INTEGER  J0ATA<12) 

EQUIVALENCE!  ATEM,  JOATA<  3  >  )  .  <  STEM .  JO  ATAC  7  > 

EQUIVALENCE! RB0TEM , JO ATA! 9  > ). ! 8ETAT , JOATA!  11)) 
EQUIVALENCE!  I  YE AR . JOATA! 4 ) ) 

EQUIVALENCE!  I  MONTH.  JOATA!  2  )  >  .  !  lOAY,  .'DATA!  3  >  ) 

C  XXX  ENO  LOCAL  VARIABLES  XXX 


C  XXX  BEGIN  CODE  SEGMENT  XXX 

C  X  t  X  CLEAR  ROW  ANO  GET  SET  TO  TYPE  REQUEST  XXX 
HAZA! 3  )>S2 
HAZA! 4 )>32 

3  CALL  PRINT! HAZA) 

HAZ!  I  )-27 
HAZ! 2  >-89 
CALL  PR  I  NT! HAZ) 

20  FORMAT!  ) 


X  X  r  CHANGE  HEAT  PULSE  DURATIONS  <  *  « 

IF  CICOM.HE  72)  goto  100 


!  I ACTVS!  1  > ,  12! 2  >  ) 

!  IHEAT!  1  ),  12!  14  >  ) 

!  IPSR.  12! 26  >  > 

!  tSR ,  12! 23  )  ) 

!  ITE.XT!  1  ,  1  ),  12!  30  >  > 
!  lACTVC!  :  ) .  12!  I SS  )  ) 
!  lAUTOC.  12!  192  > ) 


<  PPBETA!  1  >. R4!  7  )  ) 
!  PPB!  1  >. R4!  19  )  ) 

! CTEMP! 1 >. R4! 31 > ) 
C  RUALIM , RT'  43  ) ) 

<  AUTOTN - R4<  45  )  ) 

! EXINT. R4<  47 )  ) 


52 


TYPE  3Q 

39  PORMftTC  ' UPCftTE  HEAT  PULSE  DURATICM  OF  WHICH  CHANNEL  ?  ' ,  f> 

.•iCCSPT  -t'S-  I 
+0  FORMATC  I  ) 

TYPE  aa 

IF  C I . LT . 1 . OR . I . CT  S >  GOTO  3 
CALL  PRIHT<HA2A5 
CALL  PRIHTCHAZ) 

TYPE  sa 

33  FORMATC  • ENTER  NEU  HEAT  PULSE  DURATION  ',»> 

ACCEPT  Sa,tHEAT'ir> 

S3  FORMATS  13  > 

GOTO  laao 


C  «  K  «  CHANGE  FILENAME  *.  X  % 
laa  iF< icoM . HE . pa  )  goto  aaa 

TYPE  lia 

113  FQRMAT<  • UPDATE  FILENAME  OF  WHICH  CHANNEL  ?  '.*> 

ACCEPT  128- I 

laa  FaRMAT<  r  > 

TYPE  20 

IF<  I  .  LT  .  1  .  OR  .  I  .  GT  .  S  >  GOTO  5 
CALL  PRIHTCHAZA) 

CALL  PRINTCHAZ) 

TYPE  138 

133  FOR«AT< • ENTER  NEU  FILENAME  -  -  -  aAANHN  •-*> 

IP8  accept  138,  IFIL£<  1 , 1  >> IFILE<  I- 2 >-  IFILE<  I  -  3  > 

ISa  FORMATC 3A2 > 

TYPE  28 

OECOOeca,  ISa,  IFILEC  r,3).ERR-17a)lFILEA 
ISa  FORMATC 12 > 

IFC IFILEC r- 2  ). CT.  *38888 . AHO .  IFILEC I. 2 LT .  "33888 >  GOTO  1888 
1P3  CALL  PRINTCHAZA) 

call  PRINTCHAZJ 
TYPE  138 

130  FORMATC  ‘INCORRECT  ENTRY  -  TRY  AGAIN  -  AAANHN  •,»> 

GOTO  IP8 


C  XXX  CHANGE  PROSES,  TEXTS,  AND  BRIDGE  BALANCES  <  «  ;t 


288 

IF  CIC0M.NE.3e>  GOTO  488 

HAZC 1 >>2S 

HAZC  2  >-"288 

283 

CALL  PRINTCHAZJ 

TYPE  28 

TYPE  218 

218 

FORMATC 'UPDATE  PROBE  INFORMATION  OF  WHICH  CHANNEL  ? 
ACCEPT  213, I 

213 

FORMATC  r  ) 

TYPE  28 

IF  C I . LT . 1 . OR . I . CT. E)  GOTO  285 

CALL  PRINTCHAZJ 

TYPE  218,  I 

218 

FORMATC 'PROBE  IS  OH  CHANNEL  -  * 
TYPE  228, IPROBEC I > 

,  I  > 

228 

FORMATC ' PROBE  HUMBER  -  ' 

TYPE  a23,PPR88cr) 

,  13  J 

325 

FORMATC 'RSa  -  ' 

TYPE  238, PPBETAC I > 

,  Fa  2  > 

333 

FORMATC ■ BETA  -  • 

TYPE  23S,PPACIJ 

,F3 . 3> 

233 

FORMATC 'A  -  ' 

,  Fa  4  ;• 

W  f'l  tu  CJ  M  M  M  M  M  M  U  n 


’-VO 


FS  -♦  > 


Tv=E  2-»e ,  pps<  I  ;■ 

PQP.MtftTC  '  s 

TYPe  2*5  .  <  I  C3ATEC  I  ,  K  >  ,  K'- I  -  3  > 

Z-'S  FORMAT <  '  PRCEE  CP L  I  SKATED  CH  -  •  .  I  3  ^  .  I  2  ^  .  I  2  ■> 

TYPE  230- IPPRBBC I > 

230  FORMAT' ‘ 9R inCE  BALAMCE  RESISTANCE  -  -  ■.I4> 

TYPE  2S3-PPvea<I) 

333  FaRp,ftT<  .  asIDGE  BALANCE  VOLTAGE  -  -  '  ,  F3 . 2  > 

CALL  CVTTIMC  IPPTeB<  I  >  ,  KHRS  .  KM  I N  -  KSEC  -  KT  I  CKS  'J 
TYPE  3Sa -  KWHS. KMTM- KSEC 

3£0  FORMATf 'BSrOGE  SALAHCEO  AT  -  -  '-12,'  ',12-' 

TYPE  aeS- C  ITEVT^  I  - K )- K-1 . 29 > 


2S5 

FORMATC ' TEXT 
TYPE  20 

TYPE  20 

TYPE  270 

-  '.EOAEi 

2."’0 

FCRMATC 'TYPE 
TYPE  275 

•P* 

TO 

CHANCE  PROBE' > 

275 

FORMATC ' 

TYPE  280 

•  a- 

TO 

CHANGE  BALANCE' > 

2ao 

FORMATC • 

TYPE  293 

-T* 

TO 

CHANGE  TEXT' > 

2S3 

FORMATC • 

TYPE  285 

-c- 

TO 

UORK  UITH  HEU  CHANNEL'  > 

233 

283 

290 

FORMATC  '  OR 
ACCEPT  290. K 
FORMATC  A 

*■£• 

TO 

RETURN  TO  SMORGASBOA90 

IF  CK . EO . 'E' 

)  GOTO 

1000 

*  t  t  CHANGE  PROBE  HUMBER  ANO  CACtBRATrON  DATA  *  t  t 
IF  CK.HE.'P')  GOTO  330 
TYPE  295 

93  FORMA-C  'ENTER  NEW  PROBE  NUMBER  •-<> 

CALL  ASS  rCN<  3 .  ' FOB  >  PROBE . OAT ’.IS) 

READ  <3>JHUM 
tNUM-JHUMX12 

R£Aa<3><  rOATA<K:.'.K-l.  rNUM> 

CALL  CLOSe<3) 

ACCEPT  298. K 
98  F0RMAT<I3J 

J-l 

80  IFCK.ea. :OATA< J))  goto  310 

J-J412 

rF<J.LT. rHUM>  GOTO  300 
TYPE  305 

05  FORMATf'NO  PROBE  EXISTS  UITH  THAT  OES I CN AT I ON '  ) 

GOTO  217 

la  IPROSE<I>-K 

OO  328  K-1 . 12 
JOATA<  K )- rOATA<  J  > 

J-J+l 

20  CONTINUE 

PPA< I )-ATeM 
PPBC I >-BTEM 
PPR80< I )-RB0TEM 
PP9ETA< I )-BETAT 
ICOATEC I, 1 >■ I MONTH 
ICOATEs'  I.  2J-I0AY 
ICOATS< I. 35-rYEAR 
GOTO  217 

■*  t  <  CHANCE  BRIDGE  BALANCE  *  *  * 

30  IF  <r(C  ME.  '3'  >  GOTO  370 

TYPE  335 

335  FORMATC ' CURRENT  RAU  OATA 

K-1 

:F<  I  EQ .  1  1  GOTO  345 


i.. 


CURRENT  MILLIVOLTS'  > 


3-tr 


2S0 

3S3 

“sa 


3o3 

36S 

zey 

z  t 

3?0 

375 

3S0 

593 

385 

398 

C  * 
f00 

••18 

■♦28 

■♦30 


30  3-^a  j-1  -  :-i 
-fag 
OCMTISUS 
HmZAC  3 '*5 1 


HHZA<  •♦  imZZ 
TYPE  S0 

CALL  PR INT<  HAZfl  > 

:sec-Ksec 

CALL  rPOKEtf * 17«77a, «> 

JmtPSEKK ’ 1 7S772  > 

Jl-J 

rF< J . GT . *3777 ^  Jl-J-*1000a 

ppvsac  r  5-ji^3a8  3a-*7ai9 
TYPE  3Sa.  J.  PP</88<  t  > 
rORMAT<  SX,  0«^  10X>  F7  .  ■♦  > 

TYPE  355 

FORMAT<  •  PRESS  AHY  KEY  UHEM  OESTRED  BALANCE  fS  SEACHEC 
CALL  IPOKE<  ‘•■♦■♦,  "018100  .  OR  I  PEEKC  *  4-*  ;•  > 

J-rTTrNR<  ) 

CALL  IPOk:E<  "■♦■♦.  *  lE7e77  .  AMO  .  IPEEKC  "■♦4  >  > 

I  F<  J  .  CE  .  0  >  GOTO  365 
CALL  GTr«< rppTBBf r > > 

CALL  CVTTriK  IPPTSaC  I > . KHRS , KM  I H . KSEC . KTICKS  ^ 
:F<KSEC.ea. iseci  goto  3^0 

GOTO  347 
TYPE  366 

FORMATC 'ENTER  SfirOCE  RESTSTANCE  BALANCE  '.S> 

ACCEPT  367,  tPPRSSC  I  ) 

FaR«AT< 14  ) 

GOTO  ai7 

«  CHANGE  TEXT  *  *  * 

XF<K.NE.'T')  GOTO  385 
TYPE  375 

FORMAT':  ■  ENTER  HEW  TEXT  -  13345678901234567880', 

1  ' 13348678901234567890 ' > 

TYPE  330 

F0RMAT<  '  • ,  *  > 

ACCEPT  303, < ITEXT< r, K >, K-1 . 20 > 
format: 30A3 ) 

GOTO  317 

IF<K.ea.‘C*)  GOTO  205 
TYPE  390 

format: ' INVALID  ENTRY  ~  -  TRY  AGAIN’,*! 

GOTO  388 


«  CHANCE  SAMPLE  PARAMETERS  *  *  t 
if: ICON. NE. 83!  GOTO  500 
TYPE  410 

format:  • ENTER  NEW  VALUE  OF  PRESAMPLE  OURATIOH  < SEC  ! 
ACCEPT  430, IPSR 
format: 14  > 

TYPE  30 

CALL  print:  HAZA) 

CALL  print:haz:' 

TYPE  430 

format:  ■  ENTER  NEW  PRESAMPLE  FREQUENCY  :  rIZ  >  '  *  v 

ACCEPT  430- IPSP 

IP5P-60^ IPSP 

IPSR-  rpsRitsa-'  IPSP 

TYPE  30 

CALL  PR  I NT: HAZA > 

CALL  print: HA2> 

TYPE  440 
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^-3  =3RMATC  '  SNT£f?  MEU  VALUE  OF  SAMPLE  OURATION  <-  SEC  :■  ‘  .  S  : 

ACCEPT  -423  .  rSR 
TYPE  23 

CALL  PSTHTCMArAi 
CALL  PSrNT<HAE> 

TYPE  -tsa 

453  FORMATC  • ENTER  NEU  SAMPLE  FRESUENCV  C HZ >  '  . S j 

ACCEPT  423. rSP 

rsp-sa.'isp 
rsa-rsR*«a.-rsp 
care  1333 


c  * 

:*  <  CHANCE  CHANNEL  ACTI'/rTlES  «  «  T 

503 

IFC rCOM . NE . 79 )  GOTO  oOa 

TYPE  513 

510 

FORMATS ‘ ENTER  SAMPLE  CHANNEL  ACTIVITY  AS 

1  '  NUHSER  ' . «  > 

ACCEPT  Saa.IACTVS 

A  SIX  BIT  SIHAPY 

523 

FORMATC  SOI  ) 

TYPE  20 

CALL  PRINTCHAZA) 

CALL  PRINT<HAZ) 

TYPE  540 

540 

FORMATC ' ENTER  CALCULATE  CHANNEL  ACTIVITY 

I  '  SIHARY  NUMBER  '.*> 

ACCEPT  SSa-IACTVC 

AS  A  SIX  SIT' - 

5Sa 

FORMATC SO  I  ) 

CO  573  J-I.« 

IFC lACTVSC J).Ea.3>  GOTO  SCO 
lACTVSC J>-'ON' 

533 

I FC  r ACTVee  J ) . EO . 0 )  GOTO  570 
lACTVCC  J  J- • ON • 

579 

CONTINUE 

C  t  % 

t  CHANCE  ALARM  LIMITS  *  *  * 

soe 

IFC ICON . NE . 7S  )  GOTO  700 

TYPE  Eia 

Eie 

FORMATC 'TMTER  NEU  UPPER  ALARM  LIMIT  C MV > 
ACCEPT  E20.RUALIM 

323 

FORMATC F5. 2) 

TYPE  20 

CALL  PRINTCHAZA) 

CALL  PRINTCHAZJ 

TYPE  <*-30 

333 

FORMATC ’ ENTER  NEU  LOUER  ALARM  LIMIT  C MV > 
ACCEPT  S4&.RLALIM 

'  .S> 

S49 

FORMATC  Fa . 2 ) 

GOTO  1303 

C  » 

*  »  CHANCE  EXPERIMENT  REPEAT  DATA  *  *  * 

790 

IFC ICON . NE . 59  )  GOTO  033 

TYPE  710 

710 

FORMATC  •  ENTER  NEU  EXPERIMENT  INTERVAL  <■  SE 
ACCEPT  720, J 

;c  )  '  ,  *  ) 

720 

FORMATC  14  > 

S;<INT-J950 . 0 

TYPE  29 

CALL  PRINTCHAZA) 

CALL  PRINTCHAZ) 

...... '..A..,.,  ...  — 
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TYPS  ?35l 

730  FOSMOTC  '  SHTSR  N6i/  ESPEH  H-IEHT  RSPSTmONS  '  ,  t  V 

ACCSPT  7^0. ISKREP 
7-Sa  FORMAT'  r3) 

TYPE  aa 

CALL  PR INTC  HAZA  J 
CALL  PRrHT<HAZ) 

TYPE  750 

-50  FORMATC  'EMTER  NEU  TIME  UNTL  MEXT  EXPERIMENT  C  SEC  i 

ACCEPT  7S0.J 

7Sa  FORMAT','  r-f  >  I 

AaTaTN-J:«E0 
lAUTOC-a 
OOTO  1 000 


300  COMTIMUE 

:0S0  HAZ<a ^'-aaa 

RETURM 


fj  u 
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SUBROUTINE  SAMP 

THIS  SUBROUTINE  PERFORMS  THE  ACTUAL  SAMPLING  OF  THE 
DATA  AHO  CREATES  DATA  FILES  ON  DISK  FOR  EACH  ACTIVE 
CHANNEL 

C  .X  ■*  t  CONNON  VARIABLES  «  ■«  « 

LOGICAL  *!!  HAZ<  3  >- HA2A<  3  > 

INTEGER  IS<SS6> 

INTEGER*-*  I-*<12> 

REAL*-*  R-*<  S-*  > 

INTEGER  IQATA<iaBB> 

COMMON  ^HAZEL^HAZ^ HAZA 
COMMON  "VARBLE^ia^  I-*  •  R-* 

COMMON  ^TEMPS^IOATA 
C  *  *  *  END  COMMON  BLOCK  *  *  * 

C  «  *  *  BEGIN  COMMON  OESICNATIONS  «  «  * 

INTEGER  IACTVS<  6  lACTVCC  S  >.  IPROBEC  S>-  IHEATCS>.  IPPRBEto> 
INTEGER  IFILe<  6. 3  >.  ITEXT<  S,20>.  ICDATEC  S/ 3  > 

INTEGER*-*  IPPTSBCS.' 

REAL*-*  PPR80<  S  )  .  PPBETA<  6  )  .  PPA<  S  >.  PPB<  S  '  .  PPVSB<  S  > 

REAL*-*  CTEMP<S).RLTEMP<6> 

EQUIVALENCE  < I  COM , I a< ! ) > ,  < I ACTvsc I > . I ac a :  > 

EQUIVALENCE  I  I  PROBEC  !  >  .  IBC  8  )  >  ,  <  I  HEAT<  I  ' .  I  a-:  1 -f  ’  > 

EQUIVALENCE  <  I PPR88<  !  i  .  I 2< 20 > > «  C  I PSR -  I a< 2S  > > 

equivalence  < ipsp. ia<a?> > .  c isr, ia< as > > 

EQUIVALENCE  < I SP . I2< 2B  )  ) .  C  1 TEXTC 1 «  1  > .  I 

EQUIVALENCE  <  I COATE< 1 . I > /  1 2< 1 30 > > .  <  I ACTVCC  I  > ,  1 2 

EQUIVALENCE  <  I F 1  LE<  1  .  !  > ,  1 2<  1 74  )  > .  <  I AUTOC -  I a<  1 38  > 

equivalence  <  IPPTB8<  1  I4C  1  > > 

EQUIVALENCE  < PPRB0< I > . R4< I > > ,  C PPBETAC I > . R4< 7 > 1 

equivalence  < PPA< I >. R4< 13 > > .  < PPBC 1 > . R4< 1 9 > > 

equivalence  <PPVB8<  1  >.R4<  a3>'>-  <  CTEMPC  I  >  ,  R4-:  31  >  i 

EQUIVALENCE  <  RLTEMP<  1  >  ,  R4<  37  >  >  -  C  RU  AL  I  M  -  R4<  4  3  : 

C  *  *  *  END  COnON  OESICNATIONS  «  *  « 

C  *  «  *  BEGIN  LOCAL  OECLARATIONS  «  *  S 
I NTECER  I OFFST<  6 ) 

INTEGER  ITIMCCSJ 
tNTEGER»a  KTIMCa) 

INTEGER  ITIM8<6) 

INTEGER  XFILET<ia> 

INTEGER  IFILEA<7> 

INTEGER  K0ATe<3J 
INTECER*4  JTIM.ISTIM 

equivalence  <  KMON, KOATE<  I  >  >, <  KOAY- KDATEC  2  >  > 

EQUIVALENCE  <  KYEAR. KOATEC  5>-.<JTIM,KTIM<  1  >> 

C  *  *  *  ENO  LOCAL  DECLARATIONS  «  *  « 

C  «  k  *  BEGIN  CODE  SEGMENT  X  1.  x 
HAZ< 1 )>2B 
TYPE  I 

I  FCRMATC  ■ NOW  PERFORMING  SAMPLE') 

CALL  PRINT<HAZ) 

:  *  *  »  chance  last  temperatures  »  *  * 

OO  10  I>t.S 
RLTEMPC I )-CTEMP< I ) 

10  CONTINUE 

Z  t  t  t  CLEAR  ICATA  *  *  * 

CO  3  1-1  -  1900 


(0 


r 
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lOATAC I i-l 
COKTIHUe 


C  *  *  *  SET  POINTERS  INTO  lOflTft  *  %  M 
JTEMP«E+ IPSR*ISR 
K-1 

OO  ae  lal.E 

I0FP3T<  I  >-IPSR  +  a 

IP  < lACTVSC I > . EQ . a >  COTQ  88 

lOFFSTC I )-< 

K-K+JTEMP 
39  CONTINUE 

Jl-IOFFSTC  1  ) 
ja-tOFFSTC  a ) 

J3-I0FFST<  3  ) 

J-»»rOFFST<  ) 
js-taFFST<  a ) 

Jfi-IOFFST<  e ) 


C  %  ■*  t  FIGURE  OUT  SHORTEST  AND  LONGEST  HEATING  PULSES  t  *  :t 
IHERTS-iaBa 
QQ  sa  I-l>6 

IF  < IRCTVS< I >  EQ . a >  GOTO  60 
IF  < IMEATS . LE . IHEAT< I > )  GOTO  £8 
IHEATS-IHEATC I > 

Sa  CONTINUE 

INEATL«a 
OO  78  tmi,6 

IF  C  IACTVS<  I  )  .  EQ  .  a  )  GOTO  78 
IF  < IHEATL . GE. rHEAT< 1 > >  GOTO  78 
INeATL-IHEAT< I > 

?a  CONTINUE 

C  F  *  *  SET  TIMES  TO  PERFORM  PRESAMPLE  *  *  * 

ITIMP-l 

I T I ME- I PSR» I PSP+ I T I  MAT  I HEATL 
ITIMO-ITI ME- I HE ATS 
OO  7a  I-l.S 
ITIMCC I )-ITlME 
ITIMB< I J-ITIME 
79  CONTINUE 

OO  80  I-l<E 

IF  <  IACTVS<  I  }  . EQ . a >  GOTO  88 
ITIMCC  I  )-ITlME-IHEAT':  I  > 

IT1MB< I >-ITIMCC I J-IPSRSIPSP 
98  CONTIFUE 

ITIMF-ITIME+ISRSISP 

ICOUNT-1 

C  »  »  »  clear  OUTPUT  VARIABLE  *  *  * 

ipio-a  I 

C  *  »  *  SET  TIME  OFFSET  *  *  *  | 

CALL  GTIM<JTIM> 

itoff-ktim< a ) 


C  »  *  PERFORM  EXPERIMENT  *  *  * 

39  CALL  GTI«<JTIM> 

IT1M-KT1M<  a i- I TOFF 

1F< ITIM . NE. tTIMA>  GOTO  99 

ITIMA-ITIMAtI 

CALL  IPOKE<  • 176779> *  I  > 

IOATA<  J1  ?-IPEEK<  '176772  ) 

CALL  IPOKEC " 176779^ '781 > 

IOATA< ja >»IPEEK< 'l7S77a> 


L 


■  *  i*>i  jV  'iA : . 


CALL  IPOKE-'  “  17S770,  “  1001  > 

IDATA<  J3  >-IPeEK<  ■ 17S772  > 

CALL  IPOtCEC  "17«77a.  “  l-tOl  ) 
lOATAC  J-*  >»IPEE1«  "  176772  ) 

CALL  IPOKEC  *  176770.  "2001  ) 

IDATA< J3  i-IPeEK<  *  1 76772  > 

CALL  IPOKEC  “  176770.  “2-»0l  > 

IQATA<  J6  >»IPEEKC  “176772 ) 

IF< ITIM -GE . ITIME )  GOTO  220 

rF< ITIM . NE . ITI«B< 1 > )  GOTO  110 
IF< ITIM . EQ . ITIMC« 1 > )  GOTO  100 

rpro-ipro . or . *  i 

Jl»Jlfl 

ITINBC  1  )-rTIMB<  1  )+IPSP 
GOTO  110 

100  IPIO«"177776  AND  IPtO  OR.*400 

110  IF< ITIM . NE . ITIMB< 2 > >  GOTO  130 

IF< ITIM . EQ . tTIMC< 2 > >  GOTO  120 
IPIO-IPIO . OR . “2 
J2-J2+1 

ITIMBC  2  )-lTIMB<  2  >-*-IPSP 
GOTO  130 

120  IPIO«* 177773 . ANO . IPIO . OR . "1000 

130  rF< ITIM . NE . ITrMS< 3 > >  GOTO  ISO 

IF<  ITIM  EQ .  ITIMCC 3 )  )  GOTO  1T0 
tPIO»IPrO . OR. *4 
J3-J341 

ITIRBC  3  l-ITIMBC  3  )-»rPSP 
GOTO  130 

140  IP  I 0»" 177773 . AND .  IPIO . OR .  "2000 

180  IF<  ITIM . HE .  ITIM8< 4  )  >  GOTO  170 

IF<  ITlM . EQ .  ITIMC< 4 >  )  GOTO  160 
IPI0«IP10 . OR .  *  10 
1 

ITIM8<  4  >«ITIM8<4 )4IPSP 
GOTO  170 

160  rPIO-"177767, ANO. IPIO. OR. "4000 

170  IF<  ITIM . NE .  IT1MB< 3 >  )  GOTO  190 

IF< ITIM .EQ . ITIMCC 3 > >  GOTO  180 
IPIO-IPIO . OR  "20 
J3-J3+1 

ITIM8<3 )-ITIMBCS >+IPSP 
GOTO  190 

100  IPIO-" 177737 . AND . IPIO . OR . "10000 

190  IFC  ITIM  NE .  ITIMBC 6 >  )  GOTO  210 

IF<  ITIM. EQ .  ITIMCC 6 >  )  GOTO  200 
IPIO-IPIO . OR . "40 
J6-J6«l 

ITIMBCE J-ITr«SC6)*IPSP 
GOTO  210 

200  IPIO-" 177737 . AND . IPIO . OR . "20000 

210  CALL  IPOKEC  "177336.  IPIO  ) 

GOTO  90 

220  CALL  IPOKEC “ 177336 ." 100 > 

IF  C ITIM . EQ . ITIMF )  GOTO  300 

IFC  ITIM. EQ .  ITIME  )  CALL  GTIMCIST:m> 
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rcouNT-; COUNT- 1 
IFCrCOUNT.ME  GOTO  SO 

ICOUNT-ISP 
Jl-Jl-*-! 

ja-ja-n 

ja-ja+i 

JS-J3-H 
J6-J6+1 
GOTO  9a 

C  »  »  »  CLEAR  OUTPUT  LIGHTS  XXX 
300  .  CALL  IPaKE< ' 1 77336 . 0 > 

C  «  »  ;»  OUTPUT  TEMQRARY  PRESAMPLE  DATA  X  t  X 

TYPE  310 

310  FORMAT<  ■ EXPERIMENT  IS  DONE') 

TYPE  330 

3aa  FQRMAT< ' P IRST  TEN  PRESAMPLE  POINTS  OF  EACH  CHANNEL'  > 

00  393  1-1,6 

TYPE  390, <  rOATA<  J  ), J-iaFFST<  I  >,  lOFFSTf  I  >  +  9  > 

390  FaRMArC10a7> 

393  CONTINUE 

TYPE  910 

910  FORMAT<  ) 

TYPE  930, < I0FF3T< r >, I-l , 6 > 

9aa  FORMATCSiai 

C  *  *  *  OUTPUT  TEMPORARY  SAMPLE  DATA  XXX 
TYPE  910 
TYPE  330 

530  FORMATC ' FIRST  TEN  SAMPLE  POINTS  OF  EACH  CHANNEL ' > 

TYPE  910 
OO  370  1-1,6 

TYPE  360, <  IOATA<  J  >, J-IOFFST<  I  >9IPSR,  IOFFST<  I  >♦ IPSR  +  9  > 
360  FORMATCiaO?) 

370  CONTINUE 

C  XXX  OUTPUT  DATA  FILES  XXX 
IFILeA< 1 )-'F0- 
IFILEA<  a  >- •  1  ■  ' 

IFILEA<  6  >- •  .O' 

IFILSA<  7  ' AT • 

CALL  lOATEC KMON, KOAY, KYEAR > 

OO  390  I-l , 18 
IFILETC I >--30090 
390  CONTINUE 

OO  610  t-1,6 

IF  < lACTVSC I ) , EQ . 0 >  GOTO  618 
OO  600  J-1,3 
K-I*3-3-*  J 

IFILET<K>-IFILE< I.J> 

600  CONTINUE 

610  CONTINUE 

00  700  1-1,6 

IF  <IACTVS<I>  EQ . 0 >  GOTO  700 

OO  630  J-1,3 

IFILEA<  J  +  a  >-IFILE<  I . J  > 

630  continue 

CALL  ASSIGNC 3, IFILEA, 19 > 

URITE<  3  >1 ,  ISTIM, KOATE,  I .  IFILET, <  ITEXTC  I , K  > - K- I , 30 ) , 

1  IPR08E<  I > , PPRO0< I > , PPBETAC I  > , PPA<  I  > , PP9<  I  '  , 

9  f  ICOATECI ,  K  >,  <-l ,  3>. 

a  IPPRSS<  I >,  IPPTSOC  r>,PPVBB<  I  >,  IHEATC  I  >,  IPS 

3  ISR,  rsp,  <  lOATAC  K  ),  K-IOFFST<  I  >,  IGFFSTC  I  [ 


K  a. 
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reo 
c  •* 

see 

3ie 

320 


CALI-  CLOSEC  3  ) 

continue 

*  WAIT  FOR  RESPONSE  SEFORe  RETURNING  TO  WAIN  »  *  « 
IF< ICON . EQ . SS >  GOTO  820 
TVPC  800 

FORRAT< • PRESS  RETURN  TO  EXIT  BACK  TO  SNORCASBORO '  > 
ACCEPT  310,1 
FORMAT<  A  > 

RETURN 

END 


non 
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SUBROUTINE  CALC 

THIS  S'JBROUTIN  REA03  IN  THE  OES I REO  OATft  FILES 
ANO  CALLS  CRUNCH  WHICH  PERFORMS  THE 
ACTUAL  CALCULATIONS 

C  «  t  «  COMMON  variables  «  «  « 

LOCICALtl  HAZ< 3 > > HAZAC S > 

INTEGER  12036) 

INTEGER^T  :-»<I2> 

REAL»-»  R4<6T) 

INTEGER  tOATA<ia00) 

COMMON  ✓HAZEL^HAZ. HAZA 
COMMON  ^VARBLE^ 12. 14 . R4 
COMMON  XTEWPS^IOATA 
C  «  ;«  *  ENO  COMMON  BLOCK  X  *  * 

C  XXX  BEGIN  COMMON  DESIGNATIONS  *  *  * 

INTEGER  I  ACTVSC  6  ).  I ACTVCC  6  ),  I  PROBE':  6  )  .  IHEATC  S)-IPPRBSCS;' 
INTEGER  IFILE<  6.3).  I TEXTC  6 , 20 ).  ICDATEC  6.  3 ) 

EQUIVALENCE  <  I COATEC  1  .  1  >  .  1 2C  1  30  >  )  .  <  I ACTVCC  I  >  ,  I  2':  1  £3  )  > 

EQUIVALENCE  <  I F ILEC  t .  1  ) .  1 2<  1 74 > >  ,  <  I AUTOC ,  I ai  1 92  )  ) 

C  XXX  ENO  COMON  DESIGNATIONS  *  *  * 

C  X  t  X  BEGIN  LOCAL  VARIABLES  X  X  X 
INTEGER  IFILEA<7) 

C  t  :«  X  ENO  LOCAL  VARIABLES  XXX 

C  X  ;*  *  BEGIN  CODE  SEGMENT  XXX 
IFILEA<  1  )-'F0' 

IFILEA<  2  )-■  1  >  • 

IFILEA<  6  )■ '  . 0 • 

:FILEA<7)»' AT' 

HAZC 1  )-26 
CALL  PRINT<HAZ) 

DO  100  1-1.6 

IF<  lACTVCC  I  )  .  EQ  .  0  >  GOTO  100 
TYPE  3. I 

3  FORMATC 'CRUNCH  CHANNEL  NUMBER  '.I) 


DO  10  J-1.3  . 

1FILEA<  J*2  )-IFILE<  I . J  ) 

10  CONTINUE 

CALL  ASSIGN*:  3.  IFILEA.  14  ) 
REAOC 3. ERR-6a )  IDATA 
60  CALL  CL0SE<3) 

CALL  CRUNCH 

TYPE  70 

~"3  FORMAT<  ■  SACK  FROM  CRUNCH') 

100  CONTINUE 

ACCEPT  30 
FORMATC  A  ) 

RETURN 
ENO 


33 


Figure  19.  Flowchart  of  Subroutine  AUTO 
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S'JSROtJTiHE  r^UTO 

THIS  SUSROUTIHE  CftLI.  SAMPLE  AHD  CALCULATE, 
IJ-ICaEMEHTS  "ILEWAMES  AMO  ADJUSTS 
REPETITION  PARAMETERS 


■*  «  COMMON  VARIABLES  t  «  * 

L0CICAL:»1  HA2<  3  >,  HAZAC  S  > 

integer  12<a3S> 

INTEGER*-*  14<ia> 

REAL#-*  R*<6-*) 

INTEGER  lOATA<iae0> 


COMMON  y-NAZEL^HAZ- HAZA 
COMMON  ••VARBLE^'ia,  t-*,R-* 
COMMON  ^TEMPS,' lOATA 
C  *  :«  #  EMO  COMMON  BLOCK  *  *  * 


#  #  #  begin  common  designations  «  «  t 

INTEGER  I  ACTVSC  fi  >.  I  ACTVCC  S  ),  IPROBEC  6  >.  IHEATC  S  ■>  .  IPPREBC  S 
INTEGER  IP ILE<  S, 3 ).  ITEXT<  E , 20  > .  ICOATEC  S, 3  > 

C  lACTVSC  1  ).  I2<  2  >  > 

C  lACTVCC  1  >,  I2<  1S8  )  > 
c  I AUTOC ,  izr.  1 S2  >  ;■ 

C  AUTOTN  ,  R-*<  -*S  >  ■> 

<  exint,  r:*i;  *7  j 


C  «  *  S  BEGIN  LOCAL  VARIABLES  *  *  X 
INTEGER*-*  KTIME 
INTEGER  KTIN<a> 

equivalence  < KTIME, KTIM< 1  }  > 

C  *  X  *  SNO  LOCAL  VARIABLES  «  «  « 

C  X  X  X  BEGIN  CODE  SEGMENT  «  *  « 
CALL  CTIM<KT1ME> 

AUTOTU-6S336 . 0*KTIM<  1  >-»-KTIM<  2  > 
AUTOTN-EXINT 


.#  * 


equivalence 

equivalence 

equivalence 

EQUIVALENCE 

equivalence 

equivalence 

*  ENO  COMOH 


c  rcoM, tz< 1 > >. 

<  ICOATE<  1 ,  1 >,  I2C  130  >  >, 

<  IFILE<  1 .  1  >,  I2C  17*  >  > ■ 

<  1  EXREP,  I2<  133  >  > 

<  RLALIM, R*<  ** ) > . 

<  AUTOTU, R*<  *S  1  >, 
OeSICNATIQNS  XXX 


C  *  *  *  PERFORM  SAMPLE  < IF  NEEDED)  XXX 

J-0 

oo  lee  I >1,6 
IF<  IACTVS<  I  )  NE  .  0  )  J-J-»l 
100  CONTINUE 

IF<J.NE.0)  CALL  SAMP 


C  *  *  *  PERFORM  CALCULATE  < IF  NEEDED)  *  *  * 

J-0 

□O  110  I»1,G 

IF  C  I  AC7VC'  I  ?  ME  .  0  )  J-J-*l 
llO  CONTINUE 

IF  <J.NE  0)  CALL  CALC 

C  *  *  *  INCREMENT  FILENAMES  *  *  * 

OO  IBB  1-1,6 

IF< IACTVS< I > . EQ . 0 . ANO . I ACTVCC I ) . EQ . 0 >  GOTO  130 
OECOOeC  2,  120,  IFILEC  1,3) >IFILEA 
'.33  F0RMAT<I2) 

IF  ILEA- IF  ILEA-*  1 

TYPE  123, I, IFILEC I ,3), IFILEA 
-.23  FORMATC  I  ,  ■  ',A2,’  '  ,  1 3  •- 

IFC IFILEA . LT , 100 )  GOTO  163 


'nI  Ji 


«  :t  X  Aoo  lae  *  X  * 

:fileb-ifile<  i.2>-»“-»00 
TYPE  l-to-  IFILE<  I  .a>.  IPIt.E8,  IF  lUEB 
■♦0  FORPATCrtE^  ■  '.Aa,  •  '  ,  07  > 

■t  ■*  X  CHECK  GE  TO  X><X90e  «  X  « 

«  *  «  LEAVE  AT  XXX899  XXX 

IF<  IFILEB  .  GE  .  'S-ltee  >  GOTO  180 
XFILEC I<a>-IFtLEa 

0  ENCOOEC a. 170, IF1LE< I > 3 > >IFILEA 

0  FORMATC  la  > 

IF<  IF ILEA . LT.  10 )  IFILEC  I - 3>-IFILEC  I , 3  > . OR .  "60 
IF<  IF  ILEA  EQ . O )  IFILE<  I . 3  >-*30Oee 
80  COHTINUE 

X  X  X  UPQATE  AUTO  PARAMETERS  XXX 
IEXREP-IEXREP-1 
IF( lEXRCP . HE . 0 >  GOTO  200 
ICOM-- B ' 

AUTOTN-0 . 0 
I£XREP>1 

00  RETURN 

EHO 


SUeROUTIHE  DATA 

C  THIS  SUeROUTINE  ALLOWS  THE  USER  TO  EITHER  LIST 

C  AH  ENTIRE  OATA  RILE  OR  PLOT  A  FILE  ON  THE 

C  GRAPHICS  TeRMtNAL 

C  *  *  %  COMMON  VARIABLES  «  «  « 

LOCtCAL«l  HAZC HAZAC 5> 

INTEGER  I2<2SS> 

INTEGER*-*  I-*<ia) 

REAL*-*  R-*<S-*> 

I HTEGER  I O ATA<  1 Sao  > 

COMMON  ^-HAZEL^HAZ^  HAZA 
COMMON  ^VARBLEz-IS/  I-*.R-* 

COMMON  ^TEMPS^IOATA 
C  *  *  «  ENQ  COMMON  BLOCK  *  «  * 

C  t  *  «  BEGIN  LOCAL  OES I  GNAT  IONS  *  *  * 

INTEGER*-*  ISTIM 

INTEGER  KDATE<  3  >.  IFILETC  13 ),  ITEXTC  20  ) 

I NTEGER  1 CO ATE<  3  > 

REAL**  PPRBB- PP8ETA >  PP A >  PP9 *  PPVBB 
INTEGER*-*  IPPTB8 
REAL*-*  OATACIBBB) 

INTEGER  IFILEA<7) 

COMMON  ^GRIO-'IX.IY 

C  «  «  *  ENO  LOCAL  OECLARATIONS  *  *  « 

IFILEAC I }»'FO' 

1FILEA<2  >-•!•• 

IFILEA<«)«' .0’ 

IFILEA<7)-* AT' 

10  JTEMP-ITTOUR< 26 ) 

TYPE  30 

30  FORMAT<  J 

TYPE  -*0 

-*0  FORMAT<  'ENTER  FILENAME  YOU  WISH  TO  OPEN  '.»> 

ACCEPT  se.  tFILEA<  3  IFILEA<  4  >.  IFILEA<  5  > 

SO  FORMATC  A2.A2jA2  > 

CALL  ASSIGN< 3. IFILEA< 14> 

READ  <3<ERR»60>  J .  I  ST IM . KOATE .  I  -  I F ILET .  I TEXT ,  I  PROBE . 

1  PPR80.PPBETA,PPA,PP3> ICOATE. IPPRBB. IPPTBB. PPYSB- 

2  INEAT/ IPSR, IPSP> ISR> ISP, IDATA 

60  CALL  CLOSEC3> 

ICHAN>I 

6S  type  70 

70  FORMAT<  'ENTER  ‘G"  FOR  GRAPHICS  OR  'L*  FOR  LISTIHC'> 

ACCEPT  80,1 
38  FORMAT<  A  > 

IF<I.Ea.'C'>  GOTO  SOO 

C  *  *  *  LISTING  PACKAGE  «  *  * 

I-tCHAM 
TYPE  1O0,J 

100  FORMATC 'DATA  VERSION  -  -  ',!> 

CALL  CVTTIM<  I  ST  I  M- KHRS .  KM  I N.  KSEC  -  KT  I  CKS  ;■ 

TYPE  1 10, KNRS. KMIN, KSEC. KOATE 
110  FORMATC  'SAMPLE  TIME  -  -  '  .  1 2 ,  '  -  '  .  1 2 ,  '  '  ,  1 2 , 

1  '  OATE  -  -  '.12,','',:2,'-''.:z  ^ 

TYPE  1 2’3  ■  I 

120  FORMATC ' CHANNEL  OATA  RECIEVED  FROM  -  -  ' , t I 
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■ 

1  -3 

133 
ISS 
1  C3 

la3 

170 

ISO 

ISO 

230 

23S 

210 

229 
225 

230 

2-»0 

Z^S 

C  * 

390 


310 

329 

330 


TYPE  133 

P^RMPT'r  '  PCTH-'E  FILES  PT  SAMPLE  TIME'j 
TYPE  X+S, IrILET 

FOPMATC  '1-  •  .  3A2,  '  2-  ‘  .  3A2 ,  '  3-  '  -  3A2-  '  -t-  ’  ,  3A2/  '  5-  '  . 

1  3A2.  •  S-  •  - 3A2  > 

TYPE  130.: TEXT 

FORMAr< ' TEXT  -  -  ' . 20A2 > 

type  XSS.rPROOE 

FQRMAT<  • PROBE  DESIGNATION  -  -  '.ISJ 
TYPE  ISO. PPR90. PPBETA. PPA. PPB 

FORMATS ' RBO  -  -  '.F7.Z,-  BETA  -  -  '.FaS. 

1  'A  -  -  '.Fa.-*,-  a  -  -  •.Fa.'^> 

TYPE  1S5, ICOATE 

FORMAT<  ■  CALIBRATED  ON  -  •  .  12  .  ‘ .  I  2  .  ■  .  I  2  > 

CALL  CVTT  IM<  IPPT88.  XHRS,  iCMlN.  KSEC.  KT  ICKS  > 

TYPE  170. IPPRBB. PPVBO. KHRS. KMIN. KSEC 

FO««AT< ■ BRIDGE  BALANCE  .  OHMS  -  -  ’.1S.‘  VOLTS  -  -  ' . F3  2 . 

1  ■  time  -  -  • . 12. • > ‘ . ra. ‘ - 12 > 

TYPE  lao. I HEAT 

FORMAT<'HEAT  PULSE  DURATION  -  -  •-I4> 

TYPE  190, IPSR. IPSP 

FORMATC  '  PRESAMPLE  READINGS  -  -  ',1-*. 

1  • presample  PERIOD  -  -  '.iSi 

TYPE  200.  I  SR,  ISP 

PORMAT< •  SAMPLE  READINGS  -  -  '.I4. 

1  •  sample  PERIOD  -  -  ',135 

PPRflO-5 . 25^  *3777 

00  aas  K-i .  iPSR+tSR 

tF<  rOATA<K>  .  CT.  ''3777>  lOATAC  (C  >■  lOATAC  K  >- •  I  OOOO 
OATA<f<}-rOATA<  X  XPPRBO 
CONTINUE 
TYPE  210 

FOR«AT< ' PRESAMPLE  DATA  =  < M ILL IVOLTS > ' > 

00  aas  r-i.iPSR.i0 
j-r+9 

rF< IPSR-r , LT. 9 >  J-IPSR 
TYPE  a20.<0ATA<K).K-r.U> 

F0RMAT< 10F7  4 ) 

CONTINUE 
TYPE  230 

FORMATC  *  SAMPLE  DATA  •  < M I LL I VOLTS ) '  i 

DO  243  :» IPSR* 1 . IPSRTISR. 10 

J-I+9 

IF< tPSR+ISR-l . LT . 9 >  J«IPSR*ISR 
TYPE  240. <OATA<K>.K-I. J> 

FORMATC 10F7 . 4  > 

CONTINUE 
GOTO  290 

«  GRAPHICS  PACKAGE  «  X  « 

JTEMP- I TTOURC  2S  > 

■*  SET  GRAPH  SIZE 
JTEMP- ITTOURC  13  > 
rXMIN-35 
:xnAX-i0i0 
IYMIN-40 
I YMAX-740 
TYPE  310 

FORMATC • DO  YOU  WANT  DIFFERENT  SIZE  GRAPH? ' > 

ACCEPT  320. I 
FORMATC  A  > 

IFCt  NE-'Y'-,.  GOTO  343 
TYPE  330 

FORMATC  '  CONSIDER  THE  SCREEN  TO  BE  190  a’’  1 30  • 
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TYPE  93S 

37Z  FORMAT-;  'ENTER  LEFT  HORIEOMTAL  MARGIN  FROM  3  TO  99  ‘  > 

ACCEPT  S-tB  ,  I 

?-‘'3  format;  ra> 

TYPE  3-1S.  I 

S-fS  format;  '  ENTER  RIGHT  HORIEONTAL  MARGIN  FROM  '-IE-'  TO 

ACCEPT  9Sa,J 

330  format;  la) 

if; j . ot. I 3G0T0  sea 

TYPE  933 

?33  format;  •  INCORRECT  VALUES  -  -  TRY  AGAIN'  J 

GOTO  932 

3o3  I-;  IXMAX-IXMIN  >*Iy'99+IXMIN 

ixmax-int;;  ixmax- ixmin >^99  a*j>+ixMiN 

IXMIH-I 

3«a  TYPE  963 

963  format; 'ENTER  LOWER  VERTICAL  LIMIT  FROM  0  TO  99' ) 

accept  9?a.  I 
993  format; 12 i 

TYPE  973-  I 

373  format; ' ENTER  UPPER  VERTICAL  LIMIT  FROM  '-12,'  TO  99' 

ACCEPT  900.  J 

900  format; 12 > 

if; J. GT. I >caTo  99a 
TYPE  993 

983  format; ' INCORRECT  VERTICAL  LIMITS  -  -  TRY  AGAIN') 

GOTO  962 

990  I-<  IYMAX-IYMIH)RI^99-*rYMIN 

IYMAX-;  IYMAX- IYMIM  )* J^99+ IYM in 
lYMIM-I 

343  JTEMP-ITT0UR;26> 


C  t  <  *  SET  TIME  SCALE  LIMITS  Jt  S  * 

ITNEG-  tPSR*  IPSP-*-  IHEAT 
ITPaS*ISR*rSP 
ITMIHS—  ITNECx6a 
ITMAXS-ITPOS^Sa 
ISTART-  ITHEG-i-  ITMINS*6a 
ITOTAL-ITNEG-MTPOS 

TYPE  laia 

1919  format; '00  YOU  WANT  TO  EXPAND  TIME  SCALE?') 

ACCEPT  1913.1 
1913  format;  A) 

IF;l.He.'Y')  GOTO  1199 
1929  TYPE  1923. ITMIHS.  ITMAX3 

1023  format; 'ENTER  LOWER  TIME  90UH0ARY  -  -  FROM  ' , 13. '  TO 
ACCEPT  1027.  r 
1027  format; 13) 

if; r . GE. ITMAXS . OH . r . LT. ITMIHS >  GOTO  1029 

1329  TYPE  1030. I.  ITMAXS 

1330  format; ' ENTER  UPPER  TIME  BOUNDARY  -  -  FROM  ',13.'  TO 
ACCEPT  1932.0 

1032  format; 13) 

if; J.GT. ITMAXS . OR . J. LE. I )  GOTO  1089 

ITMAXS-J 

ITMINS-r 

ISTART-a 

ITOTAL-; ITMAXS-ITMINS )*Sa 
1130  JTEMP-ITTOUH; 26 ) 


99  ' 


' .  r 


•  I 


K  -t  ;  label  curve 


•«  ■«  « 
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JTSMP- trTaUR<  23  5 

rx-*a 

: V-77S 

CriLL  CPOtNT 

JTEMP-ITTCURC 13 > 

TYPE  310. ITSXT 
310  PORWpr< 3ex. 2302 > 

TYPE  320.  IFILSACSJ.  rFtLSrA<4).  IFILEACSJ.  ICHAH 
020  FORMAT<30X.  'FrUEHAMC.  ' , 3A2.  '  CHANNEL^  '.fJ 

TYPE  330.  IPHOBE.  rHEATx«0,  tHCAT- IHEAT/'EOKSa 
330  FQRMAT< SOX. • PROBE  '.12. ’  HEAT  PULSE  ',r2.‘  '.r2 

CALL  cvrrin<  rsTr?4.KHRS.XMrM.xsEc.KTrcfcs> 

TYPE  340. XHRS. KMtM. KSEC. KOATE 

340  FQRWATC sax. 'SAMPLE  TIME-  ' , 12. ' ' ' . 12. ‘  . 12 , '  DATE-  ' 

1  ra.  •>"  .  ra.  .  r2> 

JTEMP- t TTOUHC  23  5 


C  S  •;  '(  DRAU  AXIS  «  «  '« 
jTEMP-rTTOUR<  29  > 
rx-rxMiN 
rv-rvMAX 
CALL  CPatHT 
rX-tXMIM 
rV-IYMIM 
CALL  CPOIMT 
:'x»rxMAX 
lY-rvwrH 
CALL  GPOIHT 
JTSMP- r  TTOUR<  23 ) 

Z  X  t  ■*  LABEL  AXIS  X  X  « 

B»<  lYMAX-rVMIM  Jxia  .  S 
YMIO-<  lYMAX-IYMIN  >xa . 04tYMIH 
JTEMP-rTTOUR<  31 ) 

TYPE  549 

549  FORHAT<  > 

00  580  r— 3.5 
lY-rNTCYMio+ixa;-  ♦  ii 
IX-rXMIM-3S 
JTEWP- r TTOUR<  29  > 

CALL  CPOIHT 
JTEMP-ITTOUR<  31  ) 

TYPE  550. r 

550  FaRNAT<I2> 

JTE«P-rTTOOR<a9  > 
lY-rv-i 1 
rx-rx«iN-5 

CALL  GPOIMT 
IX-IXMIM 
CALL  CPOTHT 
JTEMP-tTTaOR<2a  ) 

530  COMTIMOE 

O-FLOATC  r'XMAX-rxMrfO^ITOTAL 
00  ES0  r-ITMIMS. XTMAXS 

rX-rMT<<  ISTART4<  r-ITMIHS  >X60  >XO  >-*-r.XMIH 
lY-IYMIK 

JTEMP«  rTTOOR<  29  ) 

CALL  CPOrMT 
rY-fY-3 
CALL  CPOINT 
JTSMP- ITTOURC  23 ) 

JT£MP-rTT0UR<29 ) 
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rx-rx-is 

CALL  CPOrWT 
jTEMP-rTTCURf  3l > 

Tvps  «-«a .  r 
s-*<3  rcRMAT<  ra 

JTEMP- I TTOUR<  33  > 

tjsa  ccMTiKue 

C  »  *  *  PLOT  CURVE  *  «  * 
rcouHT-a 

ITLOM- tTMtNS*Sa+  rTNEG-ISTART 
I TM  r  GH- r  TMAXS^Sa* I TWEC- r  START 
0-PL0AT<  rXMAX-rXMtM tTH TCH- I TLOU > 

B«PLOAT<  rvMAx-i  vwrN>^"iaBas 
00  Paa  r-i>tPSR-nsR 
iT-tPSPtr 

rF<  r . CT. rpsR )  iT-rTNec>rsp*<  r- tpsr  > 

1P<IT  LT . ITLOU . OR . IT  GT. ITHIGH >  GOTO  PS0 
tF<  rr  LT.  ITNEG  ANO .  IT. CT.  IPSP^IPSP  1  GOTO  733 
ix-rxMrM+rMT«< iT-iTLau>*o  > 

ICOUNTalCOUNT* 1 

IF< tCOUNT  HE. ISJGOTO  Pas 

rcooHT-a 

jTEMP-rTTOOR<  31  > 

TYPE  pas 

783  FORHATf  ) 

JTEHP- ITTOURC  38  ) 

Pas  IY-10ATA<1) 

tF< lY  GT . "3PPP >  lY-iY-'iaaae 
I Y-1HT<  YRIO  +  a*IY  > 

CALL  GPOINT 

•’aa  continue 

jterp-:ttour<3i > 

JTEHP- X  TTOUR<  a ♦ > 

HA2A<  3>>i37 
HAZAC  -i  )>7a 
CALL  PRINTCHAZAJ 


C  .«  .«  <  RAP  THINGS  UP  «  X  .< 

338  TYPE  ass 

ass  FORRAT< ■ ENTER  ‘E*  TO  EXIT  OR  -D"  FOR  PORE  DATA  ',»> 

ACCEPT  388 .  I 
388  FORNATCA) 

JTEHP-ITTaUR<  28  > 

JTE«P-ITTOUR< 31 > 

JTBNP-XTTOUR<  ES  > 

JTEMP- I TTOUR<  24  > 

JTEHP- r  rrauK< as  > 

XFCX.EQ.'E')  GOTO  318 
TYPE  383 

333  FORHAT<  ‘00  YOU  UANT  TO  REHAIN  WITH  SAME  DATA  •>  '  > 

ACCEPT  387. I 
387  FORHAT<  A  > 

IF<  I  .EQ.  -Y'  )  GOTO  S5 
GOTO  18 
318  RETURN 


laO  (0 'J  (•!  n>  Jt 
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PRCSEC  FOR 

this  PRCriRAM  MaOZFIES  THE  FILE  PROSE  OAT 

AND  uses  subroutine  ilET  TO  ACCEPT  NEU  CALIBRATIONS 

OtMeHSION  :OATAClSa0> 

COMMON  rOeSIG. IMONTH. IDAY. IYEAR, ATEM, STEM, PB0TEM- BETAT 
□IMENStON  JOATA<ia> 

EQU  I  valence  <  ATEM.  JOATAC  5  >  >.<  STEM.  JDATA<  7  •>  > 

EQUIVALENCE  RB0TEM.  JOATAC  3  >  >  ,  C  BETAT.  JDATAC  1  1  >  ) 
EQUIVALENCE  <  IDES  I G. JOATA<  1  >  >- <  I  MONTH. JOATA<  2  >  > 
equivalence  <  lOAY.  JDATA<  3  >  >.  <  IYEAR  ,  JOATA<  -»  >  5 
CALL  ASSIGNC  3.  ' FOa  PROSE . OAT'  ,  13  > 

REAO  ( 3  >  JHUM 
INUMaJNUM'Jlia 

REAO  <3)  < IOATA< K >. K»1 . INUM > 

CALL  CL0Se<3) 

TYPE  e 
FORMATC  ) 

TYPE  2 

F0RMAT<  'ENTER  “I*  TO  ADO  OR  MOOIFY  A  PROBE  V 

TYPE  3 

F0RMAT<  'ENTER  -2“  TO  LIST  ALL  PROBES') 

TYPE 

FORMAT<  'ENTER  •S"  TO  EXIT') 

ACCEPT  a. ICOM 

format<  I  ) 

IF  <IC0M.Ea.2)  GOTO  200 
IF  <tC0M.Ea.3)  GOTO. 380 
IF  < ICOM . NE . 1 >  GOTO  S 
TYPE  e 
TYPE  25 


»  «  .*  AOO  OR  MOOIFY  A  PROBE  H  *  « 

F0RMAT<  'ENTER  PROBE  NUMBER'  ) 

ACCEPT  30.IOESIG 

30  FORMAT!  13  > 

TYPE  6 
INOEX>l 

•(0  IF  C  lOESIC  .  EQ  .  IOATA<  INDEX  >  )  GOTO  150 

INDEX- INOEX.M2 
IF  < INDEX . LT . INUM >  GOTO  40 
TYPE  100 

100  FORMAT!  'NEW  PROBE  NUMBER'  ) 

TYPE  110 

110  FORMAT!  'ENTER  "Y*  TO  START  A  NEW  PROBE  WITH  THAT', 

1  ’  DESIGNATION'  > 

ACCEPT  112.  ICOM 

112  FORMAT!  A  ) 

TYPE  S 

IF  !ICOM,NE  'Y'  >  GOTO  3 
CALL  I  DATE.  I  month.  IDAY.  IYEAR  ) 

ATEM-e 
BTEM-a 
R68TEM«a 
SeTAT-0 
JNUM-JNUM*1 
INOEX-INUM+1 
INUM-INUM4.12 
CALL  GET 
OO  122  K-1  -  12 
ioata! index  )-J0ATA<  k  ) 

:hdex-imoe»  I 

C'CVTIHUE 


120 
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G3T0  5 
T'.-Pfi  13S 

FCRMftT<  'PROSE  OURKEHTL'i"  HAS  THESE  VALJEO  '  > 

OC  1S0  IC-J.13 

JOATA<  K  ;'»IOATAC  INDEX  > 

INDEX- I NOEX+ 1 
CONT INUE 
iHOEX-INOEX-ia 
'CALL  GET 
C0TQIE3 
TYPE  o 
TYPE  202 


C  *  S  r  LIST  ALL  PROBES  IN  FILE  *  t  t 

532  FOR.'<(AT<  'THIS  IS  A  COMPLETE  LISTING  OF  PPOEES  AND  DA 

DO  210  K«l.Jt(UM 

INDE:<-<»  12-11 
CIQ  233  J- 1  . 

JCATAC  J  >>iIDATA<  INDEX  > 

INDEX-  INDEX 1 

533  SCNTINUE 

Type  237. ioesig^ imcnth. ioay, iyear 

207  rORMATC  I  3 . SX .  I  2 ,  ' X •  ,  1 2 ,  ■ ^ ,  I  2 > 

210  CONTINUE 

GOTO  3 


:  t  t  X  CHECK  IF  REPLACEMENT  DESIRED  t  *  t 

330  T-.'PE  313 

310  FORMATC  'LAST  CHANCE  TO  REMAIN  WITH  EXISTING  PROBE  F 

TYPE  313 

313  FORMAT!  'ENTER  "I"  TO  REPLACE  OLD  PROBE  FILE'  > 

ACCEPT  3EO, ICCM 
320  FORMAT!  I  > 

IF  <IC0M.SE.1>  GOTO  -*00 

CALL  ASS  rCN! 3-  ' FOB ' PROBE . OAT '  ,  13  > 

'JRITEf  3  5  JNUM 
WRITSCSirOATA 
CALL  CLOSE! 3 5 
r-'PS  325 

32S  FORMAT!  'NEW  FILE  CREATED  -  -  OLD  FILE  OESTROYEO '  5 

GOTO  ■‘SB 
■NIB  TYPE  •4  10 

‘10  FORMAT'!  'NO  CHANGES  MAGE  TO  EXISTIHG  PROBE  FILE' 

300  FORMAT!  '  STILL  UNDER  CONSTRUCTION  '  > 

•420  STOP'  THAT  IS  ALL  FOR  THIS  MESS 

EHO 


CREATE  PROBE. OAT 

THIS  PROGRAM  CREATES  A'  EMPTY  FILE 

TO  as  usea  by  probec . for 

OIMENSrON  J0ATA!I2> 

CALL  ASS  rCN! 3,  ' FOB  ■ PROBE . OAT ',135 

as  ITS! 3  11 

00  5  K-I ,  12 

JOATA! K  5-B 

CONTINUE 

UR  ITS! 3  5 JOATA 

CALL  CLOSE! 31 

STOP 'FILE  CREATED' 

END 


Accept 

Month 


O  Cl  (A  A)  U  01  o  J| 
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3UBRCUTIHE  GET 

THIS  SUBRQUTIME  RLLOUS  THE  USER  TO  CHANGE  THE 
CALlBHATrOH  DATA  FOR  A  GIVEN  PRCBE 

CGMMOH  lOESIG.  IMOHTH.  ICAY,  IVEAR. ATEK>  STEM, RS0TEM. SET 
TYPE  s- roes  I G 


5 

FORWAT< 

*  PROSE 

NUMBER  - 

■  .  13-  ' 

COMMAND 

TYPE  10- 

t  MONTH 

10 

FOF.MPTC 

' MONTH 

m 

‘13.  ■ 

<  1  >  ' 

> 

TYPE  13. 

rOAY 

1  5 

FOftPATC 

'  OAY 

m 

•  .  13.  ' 

C  2  >  ' 

> 

TYPE  20- 

rVCFiR 

PCRMPT'.* 

'  year 

• 

'  .  13-  ' 

C  3  1'  ' 

> 

TYPE  23- 

ATEM 

2:*S 

FCPMATC 

'  A 

a 

•  ,  F3  H  ‘ 

C  4  >  ■  J 

TYPE  30- 

STEM 

FOPWATC 

'  B 

* 

•  , FB  4.' 

C  3  >  ■  > 

TYPE  33- 

RB0TEM 

FCRWATC 

*  RS0 

m 

•  .  FB  .  2.  ' 

<  6  ‘  > 

TYPE  40- 

BETAT 

■4C5 

F0RmAT< 
TYPE  43 

*  BETA 

m 

•  . F9  3,  ‘ 

C7  >•  > 

-♦t 

FORWAT< 

•TO  EXIT  TYPE 

<  8  >  '  > 

ACCEPT  100. rCOM 

1  00 

FORMAT< 

I  ) 

IF  < icon 
TYPE  103 

.  HE  .  1  ) 

GOTO 

113 

10S 

FORMPT< 

■  ENTER 

MONTH 

•  > 

ACCEPT  106. IMONTH 

10*3 

FORMAT< 

IB  > 

1 :0 

IF  <  rcoM 
TYPE  113 

.  NS  .  2  .■> 

GOTO 

123 

1  13 

FORMAT< 

•  ENTER 

OAV  ' 

> 

ACCEPT  lae.IDAY 

129 

IF  < ICON 

.  NE  .  3  > 

GOTO 

130 

123 

TYPE  1B3 
FORMAT< 

' ENTER 

YEAR  * 

> 

ACCEPT  106,  I  YEAR 

IF  <ICOM.HE.-t)  GOTO  l-»0 
TYPE  133 

FORMATC  'ENTER  A'  > 
ACCEPT  13S/ATEM 
FORMATC  F8  .  ) 

IP  <  ICON  .  HE  .  3  ;•  GOTO  130 


143 

TYPE  143 

F0R11AT<  'ENTER  S'  > 

ISO 

ACCEPT  136- STEM 

IF  < ICON  HE. 6)  GOTO 

160 

133 

TYPE  133 

F0RMAT<  'ENTER  R80 ' 

> 

ISO 

ACCEPT  13S.RB0TEM 

IF  <ICOM.NE.P)  GOTO 

170 

IP 

TYPE  163 

FORMAT';  '  ENTER  BETA  ' 

'  > 

1  70 

ACCEPT  136. BETAT 

IP  ^ICOM.NE.Sj  GOTO 

3 

PET'JRN 

ENO 
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